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spectral function furnishes information about the relaxation modes of the system. Along the way, 
we argue that in a typical non-supersymmetric theory with a gravity dual, there may exist a window 
of temperature and chemical potential at large N, in which a generic non-equilibrium state can be 
characterized by just a finitely few operators with low scaling dimensions, even far away from the 
hydrodynamic limit. 
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I. INTRODUCTION AND AN OUTLINE OF RESULTS 

Holography has given us a new paradigm to deal with strongly coupled systems . One of the many attractive 
features of this paradigm is that we can deal with phenomena at strong coupling in real time. 

Though there has been substantial progress in using holography to study hydrodynamics and relaxation of 

strongly coupled systems |^-^, we still lack a systematic method for studying non-equilibrium Green's functions in 
holography. The latter turn out to be extremely useful in many applications such as understanding thermalization 
^ and obtaining strongly coupled generalizations of quantum kinetic theories, to name a few. The importance of 
pursuing this direction can be readily illustrated by two examples. 

Modeling the space-time evolution of matter formed by ultra-relativistic collisions of heavy ions at RHIC and ALICE 
is a great theoretical challenge. It is equally challenging to develop reliable methods of inference for deducing this 
space-time evolution [p^ . Ultimately, it is important to not only understand how the matter thermalizes incredibly 
fast in time < 1 fm at temperature about 175 MeV (at RHIC) and subsequently undergoes hydrodynamic expansion, 
but also how hadrons and resonances are produced and transported in this so-called fireball before finally getting 
frozen chemically and thermally. Ultimately, we do infer the expansion of the fireball from the emitted hadrons. If the 
expansion of the RHIC fireball is indeed governed by strongly coupled physics, then we can expect that holography 
will not only help us in modeling the space-time evolution of the fireball, but also help us improve upon existing 
techniques like Hanbury-Brown-Twiss pion-intcrferometry used to deduce the expansion of the fireball. 

Quantum kinetic theories are already being employed to understand the dynamics of the hadron gas after the 
chemical and thermal freeze-out in the hydrodynamically expanding fireball uM- However, in order to understand 
the details of how the hadron gas comes to existence in the first place and its subsequent freeze-out, as also corre- 
lations in the emissions of hadrons, one needs quantum kinetic theories constructed using non-equilibrium Green's 
functions. Therefore, to understand such questions at strong coupling using holography, we need to develop formalism 
to systematically obtain non-equilibrium Green's functions. 

The second example pertains to holographic models of non- Fermi liquids ||l^-[l^ ^ . Holography has been successful 
in reproducing some of the features of ARPES experiments in cuprates and other strongly correlated electron systems 
- the spectral function has a pole on a momentum shell at zero frequency and also shows non-trivial scaling for low 
energy excitations. These results may be interpreted as holographic reproduction of Fermi surfaces different from that 
in Landau's Fermi liquid theory. In absence of a better way of dealing with strongly interacting fermions at finite 
density, holographic methods could provide us with useful qualitative insights. 

Nevertheless, to test such holographic models, we need to see if we can also reproduce qualitative aspects of non- 
equilibrium dynamics in strongly interacting fermionic systems. Ultimately, when the electrons are weakly interacting. 
Landau's Fermi liquid theory gives a unified way of dealing with both equilibrium and non-equilibrium phenomena. It 
is reasonable to expect that holography can do a similar job at strong coupling. Once again, we need to understand how 
to obtain quantum kinetic theory from holography, and therefore a systematic method of obtaining non-equilibrium 
Green's functions. 

There are two important issues associated with obtaining non-equilibrium Green's functions in field theory pq|. 
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^ Holographic non-equilibrium Green's functions as an aid for understanding thermalization have been studied earlier in ^ using geodesic 
approximation, etc. 

^ For interesting holographic models of Fermi liquids see |l7[. Our comments are applicable to such models as well. 
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1. There is no partition function which plays the role of generating functional of non-equilibrium Green's functions. 
As we will review briefly later, these are obtained from a generalized effective action. The effective action 
technique guarantees the full hierarchy is consistently solved and Ward identities are preserved. 

2. We cannot use conventional perturbation theory to obtain the behavior in time, like for instance, dependence 
of observables on hydrodynamic and relaxation modes. This is because usual time-dependent perturbation 
theory gives us the behavior in time in the form of a Taylor series, which fails to capture late time behavior like 
exponential decay. 

Therefore, even at weak coupling non-equilibrium field theory is hard and typically we need to make educated guesses, 
depending on the understanding of a specific system. It will be remarkable if, on the strong coupling side, holography 
can provide us with a good perturbation theory for the non-equilibrium observables we will deal with here. The lack 
of a generating functional for non-equilibrium correlation functions on the field theory side, nevertheless, makes it 
hard to use the holographic dictionary to translate such observables to the field theory side. 

The observables of primary importance are two-point correlation functions. In the vacuum, once the Euclidean 
Green's function is specified, we can analytically continue to obtain the Feynman propagator, the retarded and 
advanced Green's function etc. at equilibrium. At finite temperature too, it thus suffices to know the retarded Green's 
function, from which we can obtain other propagators like the Feynman propagator. At non-equilibrium the situation 
is different - we cannot deduce from the retarded Green's function, for instance, the Feynman propagator which will 
have independent dynamics. Nevertheless, all Green's functions can be expressed in terms of two independent, real 
observables - the spectral function and the statistical function, which we briefly review now. 

The spectral component (or spectral function) of bosonic Green's functions (in d spatial dimensions) can be defined 
as the Wigner transform (i.e. the Fourier transform in the relative coordinate r and time difference tr) of the 
commutator 



yt(w,k,x,i)= / d'^rdtre''-'^*^'-^-''^ 



'22/V 2 2 



(1) 



Similarly in case of fermionic fields, we can define the spectral component as the Wigner transform of the anti- 
commutator 



Aiio,k,^,t)= J d''rdt.e'(-*'->'-)({*(x+|,t+|),^(x-^,t-|)} 



(2) 



In both the equations above (....) denotes expectation value in a non-equilibrium state. The fermionic spectral function 



A{iij, k, X, t) = Tr 7M(a;,k,x,t) 



(3) 



The statistical function (also known as the Keldysh propagator) is defined as the Wigner transform of the anti- 
commutator of two bosonic fields 



GK(c^,k,x,i) = J d^rdtr e'('^*'-'' "-)({$(x 
or as the same of the commutator of two fermionic fields 



r tr\ f V tr 
+ 2-'+2).*(--2.'-2 



})■ 



(4) 



GK(c^,k,x,i) = / d'^rdi.e^^"*--''--) 



* 22/ V 2 2 



(5) 



All propagators can be expressed as appropriate linear combinations of the spectral and statistical functions. In 
this paper, we will be interested in the retarded correlation function in particular. It is actually more convenient to 
define the Wigner transform of the retarded correlator. In case of bosonic fields, this is defined as 



Gfl(w,k,x,t) = / d'*rd<^e*('^*--*'"")6'(t^) 



* 2 2 / V 2 2 



(6) 



Similarly for fermionic fields, the anti-commutator is used above. 

It is clear from the definitions of the spectral functions (|l|) and (||) respectively that the bosonic spectral function is 
related to the retarded correlator via A{uj, k, x, t) — — 2ImGij(w, k, x, t), while for the fermionic spectral function, the 
relation is A(a;,k, x, t) — — 2Im(Tr(7*G/j(cLi, k, x, t))). The retarded correlation function does not contain any more 
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information than the spectral function, since it is analytic in lo for a given k for every x and t. Therefore, 

Gfl(w,k,x,t = / — — (7) 

J ZTT U! — Lj' + le 

in both the bosonic and fermionic cases. 

On the other hand the Feynman propagator Gf is a linear combination of both the spectral and statistical compo- 
nents. For both bosonic and fermionic fields, prior to Wigner transform : 

Gf(x, t, y, t') = Gk(x, t, y, t') - ^ A{x, t, y, t') sign(i - t'). (8) 

Since the Feynman propagator involves the statistical function which is unrelated to the spectral function algebraically 
out of equilibrium, we cannot deduce this propagator from the retarded function in non-equilibrium states. 

At equilibrium, both the spectral and statistical functions depend only on uj and k, i.e. they are homogeneous in 
X and t, owing to translational invariance. Furthermore, they are related by fluctuation-dissipation relations : 

GK{to,k) = -i(^+nBEiuj))A{u;,k) (9) 

for the bosonic case and 

GK(t^,k) = -i(i-nFD(^))^(c^,k) (10) 

for the fermionic case, with nB^i^) — (e*^" ~ 1)^^ being the Bose-Einstein distribution and n-pD{uj) — [e^^ -f 1)^^ 
being the Fermi-Dirac distribution. 

Away from equilibriuin, the statistical and spectral functions follow a coupled set of equations which were first found 
by Kadanoff and Bayni These equations are not so easily tractable in field theory even at weak-coupling, however 
educated guesses lead us to standard kinetic equations like the Boltzmann equation with quantum corrections. We 
will skip issues involving renormalization etc. and simply mention here that they can be dealt with efficiently at the 
level of the effective action. 

The spectral function, especially for fermions, is directly measurable by ARPES like experiments. Usually it is the 
equilibrium spectral functions that are measured experimentally, so that we need be concerned with their dependence 
on frequency and momentum only. Recently however, there have been time-resolved ARPES experiments in which non- 
equilibrium time-dependent spectral functions have been measured in approximately spatially homogeneous situations 
and their dependence on frequency, momentum as well as time have been obtained (see, for example, time-resolved 
ARPES across the metal-insulator transition in jl9|). Conceptually, when integrated over frequency at a given 
momentum and at a given point in space-time, the spectral function gives the space-time dependent density of 
states. The spectral function thus reveals the non-equilibrium structure of the effective phase-space of quasi-particles 
(provided we do have well defined quasi-particles). 

The statistical function, on the other hand, carries complementary information about how quasi-particles (whenever 
they can be defined) are distributed in phase-space and time and can be indirectly measured. For instance, in the 
case of a single species of fermions, the conserved current is 

j^(x,t) j dw(i''fcTr(^7''GK;(w,k,x,t)) + constant, (11) 

where q is the conserved charge of the fermionic field, and the constant is independent of the state and required to 
provide an infinite subtraction which produces a finite result. In the so called quasi-particle approximation, we can 
assume that the statistical function is peaked only when w is on-shell, so that it reduces to the standard phase-space 
distribution which follows the semi-classical Boltzmann equation in certain limits. 

This completes our very brief review of the spectral and statistical functions respectively. In this paper, we would 
like to obtain the non-equilibrium retarded function holographically. Our focus will be on the retarded function 
because we can compute it using linear response theory even in a non-equilibrium state. The holographic dictionary 
enables defining the source and expectation value of an operator in any arbitrary state. Therefore, we can avoid issues 
associated with the lack of a generating functional for non-equilibrium correlation functions. 

To be specific, we would like to achieve the following : 



1. to evaluate the retarded correlation function and the spectral function in non-equilibrium states. 
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2. to find space-time dependent shifts in the energy and spins of quasi-particles in the non-equihbrium medium, 
and 

3. to obtain the space-time dependent shifts in energy per particle and spin orientation at the holographic Fermi 
surface. 

With respect to the last point, we will reproduce a strongly coupled version of what is expected from Landau's 
Fermi liquid theory, as reviewed later. The second objective is justified on the grounds that it is known that in 
non-equilibrium states, the effective masses of quasi-particles become space-time dependent (via an inhomogeneous 
temperature distribution for instance, or an inhomogeneous distribution of the velocity field as discussed later) . We 
will succeed in all these objectives for scalar and fermionic operators. 

This paper thus finishes only half of the complete formalism required to obtain quantum kinetic theory from holog- 
raphy. We do not address the information contained in the statistical function and how to obtain it holographically. 
Work in the latter direction will appear in . These issues will be complicated by the fact that we are dealing with 
composite operators in holography and we leave this for future study. We note here that there has been previous 
work where the equilibrium statistical function has been defined holographically in a consistent manner , based on 
the correspondence between the generating functional of field-theoretic correlation functions and a suitable partition 
function of quantum gravity. However, these cannot be readily generalized to non-equilibrium states because of the 
lack of a generating functional for non-equilibrium correlation functions as observed before. 

The key result in this paper will be the development of perturbation theory of scalar and fermionic fields in 
holographic duals of non-equilibrium backgrounds. At equilibrium, the incoming boundary condition mimics causal 
response in field theory and suffices to define a well-defined linear response theory holographically However, 
the incoming wave boundary condition does not suffice to give well defined linear response theory in non-equilibrium 
states. This can be briefly demonstrated as follows. 

Suppose we have a non-equilibrium background in which a hydrodynamic mode with momentum k(i,-| has been 

excited. Let the source of the operator at equilibrium be j'^^^ (x, t) and the expectation value be O'^^^ (x, t) which can 
be read-off from the profile of the field (r, x, t) in the bulk. The non-equilibrium bulk contribution can be denoted 
as 3>^'^^(k(h),r, X, i) and this gives contribution to both the source J^^^ and expectation value O^^^ of the operator. 
The full retarded function can be obtained from : 

0(")(x,i) -HO(iHk(h),^x,i) ^ ^ 

JW(y,t') + J'^^>(k(h),r,y,t') 

where C is a constant which depends on the action and has been set to unity here. However, the general solution for 
^^^^ will have : 

i) two homogeneous solutions which are incoming and outgoing at the horizon respectively and, 

ii) a particular solution which will be completely determined by the hydrodynamic background perturbation and the 
equilibrium solution 

This particular solution will contribute to both O*^^) and J*-^'', as will the homogeneous solutions. The incoming 
boundary condition will set the coefficient of the outgoing homogeneous solution to zero. The coefficient of the 
homogeneous incoming wave solution is left arbitrary. At equilibrium, this arbitrary coefficient cancels between the 
numerator and denominator, but at non-equilibrium we have an extra coefficient from (f^^-* and therefore is 
ill-defined. 

In this paper, we show that careful treatment of regularity of the solution at the horizon implies that the coefficient 
of the homogenous incoming solution should also be zero in presence of background quasinormal modes. This will 
allow us to put forth a well-defined prescription for obtaining the non-equilibrium retarded Green's function and 
spectral function holographically. In fact, the prescription can be precisely stated in a manner which is independent 
of the non-equilibrium state. Thus, holography gives a very well-defined perturbation expansion of non-equilibrium 
observables which can be understood in an universal manner. 

The organization of the paper is as follows. In section II, we give a general review of holographic duals of non- 
equilibrium states. Though most of this section is a review, the explicit metrics for charged hydrodynamics and 
homogeneous relaxation in section II. D in ^^5*4 are new as far as we are aware of the literature. The key point in the 
discussion in section II. B however, to the best of our knowledge, is novel. Here we argue that in a non-supersymmetric 
theory with a gravity dual, there may exist a window of temperature and chemical potential at large N, in which a 
generic non-equilibrium state can be characterized by just a finitely few operators with low scaling dimensions even 
far away from the hydrodynamic limit. We also point out that there are surprising similarities with solutions of the 
Boltzmann equation on the weak coupling side, which we review in section II.A. 

In section III, we develop the formalism for obtaining non-equilibrium retarded Green's function and spectral 
function holographically in the approximation where the background fluctuation is linearized i.e. when the non- 
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equilibrium state is studied in the linearized approximation. An interesting result is that we can read off the relaxation 
modes in the background by measuring the non-equilibrium spectral function. 

In section IV, we compare our holographic approach with field theory. We also make a comparison with Landau's 
Fermi liquid theory regarding non-equilibrium dynamics of the Fermi surface. Furthermore, we obtain a holographic 
prescription to calculate space-time dependent non-equilibrium shifts in the energy and spin of the quasi-particles. 

In section V, we show that our prescription for the holographic retarded Green's function rciadily generalizes when 
we take non-linearities in the dynamics of the variables characterizing the non-equilibrium state into account. 

Finally, in section VI, we conclude by pointing out interesting issues that could be addressed numerically. 

II. ON NON-EQUILIBRIUM STATES, THEIR HOLOGRAPHIC DUALS AND QUASI-NORMAL 

MODES 

An equilibrium state can always be characterized by a few macroscopic variables related by an equation of state. 
The distribution functions of particles, density of states, expectation values of operators, Green's functions, etc. 
depend on these macroscopic variables. We also know, in principle, how to calculate the equation of state relating the 
macroscopic variables of equilibrium states. Most importantly, we know in principle how to calculate the dependence 
of the observables in the underlying field theory on these variables characterizing equilibrium states. 

The most pressing problem in dealing with non-equilibrium states is that, typically even at the coarse-grained level, 
we need an infinite number of macroscopic variables to characterize them. These variables also depend on space and 
time. Aside from taking recourse to a kinetic approximation, which is typically uncontrolled (but intuitively well- 
motivated) from the point of view of the exact field theory, we usually do not know how to obtain the equations of 
motion of these macroscopic variables (thereby generalizing the notion of equation of state applicable at equilibrium). 
It is also not clear how to relate observables in the field theory to the macroscopic coarse-grained non-equilibrium 
variables. 

Here, we will address thesc^ issuers from the point of view of holography. Firstly, we will identify a special sector 
of non-equilibrium states which can be described in terms of a finite number of operators of low scaling dimensions 
in kinetic theories. These states exist for any value of the coupling at least in the kinetic approximation. Then we 
will argue holographically that these states also exist in the exact field theory and are generic at strong coupling and 
large N after a microscopic time-scale, irrespective of the initial condition. We will further discuss how solutions in 
gravity describe such non-equilibrium states. 

A. Conservative states in the kinetic approximation 

Let us first look at the kinetic approximation in some details. In particular let us analyze the Boltzmann limit 
which is valid typically when, nlf^^j^ is small, where n is the typical number density, Imfp is the mean free path and 
d is the number of spatial dimensions. 

Boltzmann equation describes the dynamics of particle-distributions in phase space. It can be reduced to local 
dynamics of the infinite number of moments of the phase-space distribution of particles f^^^ (x, p, t) of a given species 
s. These moments are 

41....„(x,i) = P,.P,.....p,J^'\^,^,t) (13) 

where p'^ is the d -\- 1-momentum with being on-shell energy for each species s. 
A conserved current (for instance the baryon number current) is given by : 

j,{^,t) = Y^qJ^ p^/W(x,p,0, (14) 

where qs is the charge (for instance baryon charge) of the s— th species. 
The energy-momentum tensor is given by 

V(x,t) = E/ ^ p,p.f^^K^,i?,t). (15) 
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Thus we see that the energy-momentum tensor and conserved currents are parametrized by a weighted sum of first 
few moments of the quasi-particle distribution functions. 
Three comments are in order here : 

1. The Boltzmann equation has no dependence on temperature or non-equilibrium parameters. The latter 
parametrize the solutions. The thermal Bose-Einstein or Fermi-Dirac distributions are exact solutions of the 
Boltzmann equation. In absence of external fields, Boltzmann's H-theorem indicates all solutions finally equili- 
brate into thermal Bose-Einstein or Fermi-Dirac distribution. 



2. The integrals involved in collision terms on the right hand side of the Boltzmann equation (see eq. (106) for 
weakly interacting electrons) have divergences coming from phase-space volume. To regulate these divergences 
one can put a IR-cutoff corrsponding to the thermal mass of the quarks and gluons with temperature being the 
final equilibrium temperature |24|. The dispersion relations are also accordingly modified. 

3. In the dilute limit the effect of the interactions is taken into account via an effective thermal mass. Thus the 
energy-momentum tensor takes a free particle form with an effective thermal mass. 

It can be shown that the higher velocity moments parametrize the flow of the flow, the flow of the flow of the flow, 
etc. of charge, energy and momentum. For example, if we define : 

Pt^PvPpf^''' (x, p, t), (16) 

then the heat-current is Sf^ = Sf^^pi]'^'' . 

The Boltzmann equation can have solutions where the partial conserved currents are are all proportional to 
each other. This happens precisely when chemical equilibrium is achieved, and in fact any arbitrary solution achieves 
chemical equilibrium after sufficiently long time. In that case, we can define a four- velocity field u'^ and charge density 
p such that : 



The energy-density e is 



e ^ tp^u^'u\ (18) 



The hydrodynamic variables are e, p and . We can define temperature T and chemical potential /i fields in terms 
of e and p by using the equation of state of the full system at thermal and chemical equilibrium locally. 

There are special solutions of the full non-linear Boltzmann equation, known as normal solutions in the literature, 
which are purely hydrodynamic These solutions are such that all the moments flf....fj,^ of the phase-space 

quasi-particle distributions of various species are algebraic functions of just the hydrodynamic variables it^, T and 
p, and their spatial derivatives in the local inertial frame co-moving with u^. The full phase-space distributions can 
thus be characterized uniquely by the hydrodynamic variables. Furthermore, any arbitrary solution of the Boltzmann 
equation can be approximated by an appropriate normal solution after a sufficiently long time. 

The hydrodynamic equations can be derived from the Boltzmann equation; these are the Navier-Stokes equation, 
charge diffusion equation and Fourier's law of energy transport with systematic higher derivative corrections. The 
shear viscosity, charge diffusion constant, thermal conductivity and all the higher order transport coefficients can be 
obtained from the relevant Boltzmann equation specified by the dominant collision processes. 

These solutions can be further generalized to what were named conservative solutions In such solutions, the 
various moments fjfi....ii„ are algebraic functionals of p, (or equivalently the conserved current j^) and the energy- 
momentum tensor tf^u, and their spatial derivatives in a local inertial frame co-moving with . Thus the full solution 
can be specified by and j'^j. In such solutions the energy- momentum tensor is not necessarily hydrodynamic. 
Furthermore, any solution of the Boltzmann equation reduces to an appropriate conservative solution after sufficiently 
long time, and the latter reduces to an appropriate normal solution after the relaxational time scale. The first claim 
follows from the fact that the independent dynamical parts of higher moments of the quasi-particle distributions decay 
faster compared to the non- hydrodynamic relaxational mode of the energy-momentum tensor [ p6[ . 

The energy-momentum tensor tf^i, and the conserved current (or equivalently the charge density p and the velocity 
u^) follow a closed system of equations in conservative solutions of the Boltzmann equation. This gives a systematic 
generalization of phenomenology beyond hydrodynamics to include processes like relaxation. These phenomenological 
equations have been obtained in . 
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Obviously, the existence of normal and conservative solutions of the Boltzmann equation can be seen at the linearized 
level and provides a method to obtain good approximations to the transport coefficients and relaxation parameters. 

Thus, in the semi- classical kinetic limit captured by the Boltzmann equation, an arbitrary non- equilibrium state can 
be approximated by a conservative state whose dynamics is given by the conserved current and the energy-momentum 
tensor even away from the hydrodynamic limit. This approximation is reliable after a microscopic time-scale which is 
shorter than the leading non-hydrodynamic relaxation mode, i.e. the time scale of local thermalization. 

The quasi-particle distribution is said to have locally thermalized when it can be characterized well by space-time 
dependent parameters of equilibrium distribution. Afterwards, hydrodynamics takes over and the system equilibrates 
globally. In a generic solution of the Boltzmann equation, we thus have three time scales. The first time-scale is the 
time for chemical equilibration tchem after which inelastic collisions effectively cease, the second time scale is t^ons 
after which an approximation by an appropriate conservative solution becomes valid, and the third time scale is after 
which the hydrodynamic approximation is valid and is also the time scale of thermalization ttuerm ■ The hierarchy is 

^chern ^ ^cons ^ ^therm- 

The conservative solutions of Boltzmann equation describe the dynamics of both thermalization and hydrodynamics 
in an unified framework in the Boltzmann limit. 

We note that there is no scale which parametrically separates the dynamics of the non-hydrodynamic part of the 
energy-momentum tensor and conserved currents from that of other relaxation modes. Thus we may argue that 
even if conservative states exist beyond the Boltzmann limit, they may not be typical non-equilibrium states after 
microscopic times as in the Boltzmann equation. The typicality is just a special feature of the Boltzmann limit. 

In fact, once we go away from the dilute limit necessary for the Boltzmann equation to be reliable or consider genuine 
quantum dynamics (not just quantum statistics), the typicality of conservative states will no longer be preserved. The 
conserved currents and energy-momentum tensor do not seem to capture generic dynamics beyond the hydrodynamic 
limit. Conservative solutions may exist beyond the Boltzmann approximation, but only in the purely hydrodynamic 
limit can they approximate a generic state. 

We will argue that if a theory has a holographic dual, then in certain phases in the large N limit, the dynamics 
can indeed be captured by just the conserved current and energy-momentum tensor generically, after a microscopic 
time-scale which is much shorter than the time-scale for local thermalization. In such cases, the conservative state 
can indeed capture generic non-equilibrium dynamics even far away from the hydrodynamic limit. 

B. Holographic duals of non-equilibrium states and typicality at strong-coupling 

Holography maps a field theory to a quantum theory of gravity in one extra spatial dimension. It further states 
that in the large iV and strongly coupled limit, the dual theory of gravity reduces to a classical theory. Therefore, 
in this limit states of the field theory are dual to solutions of the classical theory of gravity which are regular in an 
appropriate sense. Furthermore, every operator is dual to a field and the expectation value of an operator in a state 
can be obtained from the asymptotic behavior of the dual field in the corresponding gravity solution. 

The question of which operators matter in characterizing states in the large N and strong coupling limit can be 
seen from the masses of the dual fields. The mass of the field is related to the scaling dimension of the dual operator. 

The large N limit in the {D dimensional) field theory side is the limit when the scale I, corresponding to asymptotic 
curvature radius of the {D -\- 1 dimensional) space-time, is large compared to the effective Planck scale Ip (in _D + 1 
dimensions) on the quantum gravity (string theory) side of the holographic correspondence. The strong coupling limit 
on the field-theory side is the limit when the length of the fundamental string Is is small compared to the asymptotic 
curvature radius I on the quantum gravity side. The first condition Ip « I allows us to consider the classical limit of 
gravity. The second condition Is << I allows us to ignore the massive stringy fields corresponding to higher excitations 
of the fundamental string. 

Nevertheless, string theory is a theory in 10 dimensions. So, there has to be a compact space of 9 — D dimensions 
on top of the D -\- 1 dimensional non-compact coordinates. The condiitons Is << I and Ip << I, i.e. strong coupling 
and large N limit in the field theory side allows us to decouple the massive stringy modes whose masses scale like l~^ 
when Is and Ip are small compared to I. Thus from the ten-dimensional viewpoint we are left with just the massless 
fields which include the graviton and gauge fields. However, the compactification over the compact 9 — D dimensions 
still creates a tower of Kaluza-Klein fields which are dual to operators with possibly small scaling dimensions if the 
typical size of the compact dimensions is of the same order as the asymptotic curvature radius I . 

In a supersymmetric set-up , the typical radius of the 9 — 1? dimensional compact space is indeed of the same 
order as the D -\- 1 dimensional asymptotic curvature radius /. Therefore, in the strong coupling and large N field- 
theoretic limit, the Kaluza-Klein spectrum still plays a role in characterizing states. In fact, these Kaluza-Klein fields 
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are dual to chiral primary operators and their descendants. Therefore, a prediction of the holographic correspondence 
is that at large N the scaling dimensions of the chiral primary operators do not deviate much from the weak coupling 
limit. 

Despite the presence of the Kaluza-Klein spectrum, it is known that almost all known supergravity theories in 
10 dimensions admit consistent truncation at the classical level to gauged supergravity in D + 1 dimensions when 
dimensionally reduced over the appropriate 9 — 1? dimensional compact space. The D + 1 dimensional graviton is 
dual to the energy-momentum operator on the field-theory side and the D + 1 dimensional gauge fields are dual to 
the conserved currents with the global symmetry groups being gauged in the gravity side. 

One can also show that all solutions of D -I- 1 dimensional gauged supergravities which thermalize to black branes 
with regular future horizons can be characterized uniquely by the expectation values of the energy-momentum tensor 
and conserved currents of the dual states These solutions thus correspond to special non-equilibrium states - 
namely the strongly coupled version of the conservative states which can be characterized by the energy-momentum 
tensor and conserved currents alone. The parameters of phenomenological equations for the energy-momentum tensor 
and conserved currents which generalize hydrodynamics should now be obtained from gravity and not from kinetic 
theories valid at weak coupling 1^-^. Evidence that the solutions of pure gravity in particular, which have regular 
future horizons, can be interpreted as conservative states has been presented in j?] for the special case of homogeneous 
relaxation. It has been proved that regularity at the horizon gives an equation of motion for the non-hydrodynamic 
energy-momentum tensor with precise coefficients for this case. 

Furthermore, such conservative states should also exist holographically away from the strong coupling and large 
N limit, since the dual solutions in gravity can be constructed by perturbatively correcting the gauged supergravity 
solutions in I'^/P and 1/7V^. Nevertheless, in the known supersymmetric cases these solutions are always special and 
not typical even in the strong coupling and large N limit, because the intrinsic dynamics of Kaluza-Klein modes are 
absent in these solutions. 

The situation can be expected to be very different in non-supersymmetric cases. There is no analogue of chiral 
primary operators and typically we do not expect that quantum corrections to scaling dimensions of operators will 
be small at strong coupling, unless these are suppressed because of symmetries. 

In order to use our intuition obtained from well studied examples with the field theory being conformal and 
supersymmetric, we will need to focus only on a certain window of temperatures and chemical potentials, such that : 

1. the effective coupling is strong, 

2. the beta function is vanishing or approximately so, i.e. the system is close to a critical point, and 

3. there are no new emergent symmetries at the critical point other than the (exact or approximate) full conformal 
symmetry. 

Furthermore, we also require that the large N approximation is valid, or useful for qualitative understanding. Probably, 
all these requirements could be satisfied for the fireball at RHIC near temperatures of 175 MeV and small baryon 



charge densities as supported by lattice data [29 . We can also hope that the strange metallic phase of strongly 
correlated electron systems will satisfy these requirements in a window of temperatures and chemical potentials. 

We note that certain examples of non-supersymmetric holography have been proposed in the literature |30[ . How- 
ever, in these special examples, infinite number of gauge symmetries appear in the bulk at large N, implying infinite 
number of global symmetries in the dual field theory. Our observations below will not be necessarily true in such 
cases ^. 

In case of a typical non-supersymmetric theory with a gravity dual, at temperatures and chemical potentials such 
that the system is close to a strongly coupled critical point, we expect there will be a few operators whose scaling 
dimensions will be small. We observe that the scaling dimensions depend on the scale through the coupling and hence 
also on the phase of the theory being considered which is parametrized by the temperature and chemical potential. 
The relevant operators with small scaling dimensions in the window of temperature and chemical potentials considered 
here can be expected to be 

1. the energy-momentum tensor, 

2. the conserved currents, and 



^ Despite these not being Cauchy data from the gravity point of view, this holds if the geometry corresponds to regular perturbations of 
a black brane at late time |28| . We also note that the consistent truncation to pure gravity does not involve separation of scales. This 
simply reflects the fact that the conservative states are not typical states in these examples. 

* The examples [B0|| are also not stringy and so far well defined only in the large N limit, i.e. only when the theory of gravity is classical. 
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3. order parameters of spontaneous symmetry breaking. 

Therefore, the operators dual to the Kaluza-Klein modes of gravity are expected to have large scaling dimensions very 
simlar to those dual to the stringy modes. If this expectation is true, the typical scale of the compact dimensions 
should be of the same order as Is and not I. 

For instance, in the case of QCD, the relevant operators with small scaling dimensions in the conditions of RHIC 
can be expected to be 

1. energy-momentum tensor, 

2. the baryon number current, 

3. the approximately conserved SU{3)l x SU{3)jf_ flavor symmetry of the light quarks, and 

4. the order parameter of chiral symmetry breaking having zero baryon number, transforming as (3l,3_r) under 
the flavor symmetry group and with scaling dimension approximately 3. 

The dual massless fields on the gravity side should be 

1. the graviton, 

2. a U(l) abelian gauge field, 

3. SU{3)l X SU{3)ii non- Abelian gauge fields, and 

4. a neutral scalar field transforming in the (3l, Sr) representation of the non- Abelian gauge group and with mass 
approximately given by m? = —3/P ^. 

Such a holographic model for QCD has already been proposed in However, our arguments above show that such 
models can be considered more seriously in the conditions of RHIC. In fact, for RHIC conditions we also do not need 
the hardwall cut-off proposed in these models to achieve confinement, as the mass gap is expected to become very 
mild at temperatures close to 175 MeV and for small baryon number densities. 

Furthermore, if the temperature is higher than 125 MeV, chiral symmetry is expected to be restored, so that the 
profile of the bulk scalar field dual to the chiral symmetry breaking order parameter will be stabilized by a potential. 
Therefore, only the conserved currents and energy-momentum tensor can characterize non-equilibrium dynamics at 
large N and large 't Hooft coupling A for temj^eratures above 125 MeV. The other fields in the holographic dual 
should have masses which grow like 1/1^ i.e. l/A^ , and thus are expected to be effectively decoupled from the classical 
theory. 

The correlation functions of the non- Abelian gauge fields in the gravity backgrounds which thermalize to a black 
brane are all we need to construct quantum kinetic theories of production and freeze-out of axial and vector mesons 
(and resonances) in the expanding fireball holographically. The interpretation of poles of correlation functions of these 
gauge fields in terms of mesons has been given in |3^. Using the methods to be described later, we can obtain the 
non-equilibrium corrections to these mesonic poles systematically. 

Let us estimate the relevant time scale at strong coupling after which the conservative solutions become relevant. 
This in the dual gravity description is given by the mass of the lightest stringy field or Kaluza-Klein mode. According 
to the discussion above, the time scale should be 0{X^i) in a non-susy conformal theory at strong coupling. After 
such a time-scale, we may expect that the massive fields in gravity will decay and the relevant dynamics will be 
described by the metric, gauge fields and the light fields dual to order parameters of symmetry breaking relevant at 
the critical point. Thus decay of a massive field in gravity can be interpreted as transition to a conservative state 
at strong coupling where the dynamics is governed by the energy-momentum tensor, conserved currents and order 
parameters alone. 

We conclude in a typical non- super symmetric theory which has a holographic dual, in a window of temperature and 
chemical potentials such that the dynamics is strongly coupled and approximately conformal, all non- equilibrium states 
can be characterized by just the energy-momentum tensor and conserved currents (and order parameters of spontaneous 
symmetry breaking if any), irrespective of the initial conditions, after a microscopic time-scale which scales with the 
coupling A like l/A^ in the large N limit. In other words, conservative states are typical states irrespective of the 



^ As the chiral symmetry breaking order parameter is (q^q^), it has approximate mass dimension of 3. Moreover, QCD being asymptotically 
free, the dual boundary condition will be approximately AdS^-lihe as well. Then we can use the standard relation for AdS^ for mass of 
the field m and the scaling dimension of the dual operator A which gives = —3/1^ when A = 3. 
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initial conditions after a microscopic time-scale much smaller than the time-scale of thermalization in the strongly 
coupled and nearly conformal phase at large N. 

If the above arguments are indeed relevant for QCD and strange metals in a window of temperature and chemical 
potentials, we have a unique opportunity to understand non-equilibrium dynamics with only a finitely few operators 
in this special phase of these theories. As conservative states will be typical non-equilibrium states, we can use general 
phenomenological equations for non-equilibrium dynamics as proposed in , and also hope to construct a general 
theory of kinetics and fluctuations to connect to experiments as we want to do here and more completely in the future. 

If the above arguments fail, the reasons should certainly be deep. In that case, we also need to know how to 
generalize non-equilibrium holography beyond the sector of conservative states sufficiently so that we can describe a 
typical non-equilbrium state. 



C. Quasinormal modes 

The thermal states in the field theory at large N and strong coupling are captured by black brane solutions of clas- 
sical gravity holographically. In the linearized limit, the non-equilibrium fluctuations are captured by the linearized 
equations of motion of gauge field and the metric fluctuations about the black brane background. These fluctuations 
are dual to perturbations of the energy-momentum tensor and conserved currents about thermal equilibrium. Fur- 
thermore, these fluctuations should satisfy the incoming boundary condition at the horizon and Dirichlct boundary 
condition asymptotically Thus they are quasinormal modes capturing intrinsic fluctuations in the dual field 

theory which can exist in absence of sources and provide good approximation to a typical non-equilibrium state close 
to equilibrium at strong coupling and large N. 

There is, however, a significant difference between the linearized Boltzmann limit and the quasinormal mode 
approximation of solutions of gravity. Instead of a finitely few decay modes on top of the hydrodynamic mode, we 
have an infinite tower of quasi-normal modes. The reason that we do not have an infinite tower of modes for the 
energy-momentum tensor perturbations in the Boltzmann equation is that it has only one time derivative (which in 
a Lorentz- invariant language is the derivative along the local velocity field). Quantum corrections to the Boltzmann 
equation are known to result in an infinite number of time derivatives, and it is not hard to see this will produce an 
infinite number of decay modes as well. 

Wc will now obtain the phenomenological form of the non-equilibrium energy-momentum tensor and conserved 
cuurent. Instead of stating in a Lorentz-invariant way, we will state the form of the energy- momentum tensor in 
the frame where the dual thermal state is at rest, i.e. the laboratory frame. It is convenient to define the velocity 
perturbation du{x, t) such that the velocity field is co-moving with the energy-flow, instead of the charge-flow as done 
usually in the Boltzmann limit. Thus the non-equilibrium energy-momentum tensor thus takes the Landau-Lifshitz 
form in the global co-moving frame : 

too = e{T,^i) + — ^^^dT(x,t)H — — ^^l(:>c,t), 

= p{T,iJi)5,, + (M£llA5T{x^t) + ^P|liilj^(x,t))% +^,,(x,t). (19) 

Above p is the pressure and Hij is the shear-stress tensor. The shear-stress tensor can thus be defined as the dissipative 
part of the energy- momentum tensor or the spatial components of the energy-momentum tensor not in local equilibrium 
in the co-moving frame. The conserved current takes the form : 

Jo = piT,n)-\ — dT(x,t)H — — (5/4x,t) + i^o(x,t), 

J, = p{T,^i)Su^{x,t) + ly,{x,t). (20) 

Above Ui is the dissipative part of the consevred current or the spatial components of the current away from local 
equilibrium in the co-moving frame. However, as the co-moving frame is aligned with the energy fiow, the charge can 
have a non-equilibrium part by itself. This is vq- 

In order to have conformal invariance, we should further have 

e{T, fi) ^ d p{T, fx), 6e = d Sp, tt^jiJij = 0, (21) 
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with d being the number of spatial dimensions in the field theory. Above Se and Sp denote change in energy density 
and pressure due to change in temperature and chemical potential. From now onwards, we will be interested in the 
specific case when the field theory is conformal, so that on the gravity side we will be using asymptotically AdS 
boundary conditions. 

The shear-stress tensor and the dissipative part of the current can be split into hydrodynamic parts tt^^^ and v^'^'' 

respectively which are functions of the hydrodynamic fields ST and du, and non-hydrodynamic parts tt^^^"^ and t'f"'^'' 
respectively which cannot be parametrized by hydrodynamic variables alone. On the other hand, i^q does not have 
any purely hydrodynamic part. 

In the case of a conformal field theory, at the linearized level, 

(h) , (nh) (h) , (nh) 



= -T^{T,ii)\^ — — — d^5T+ — — — d,5n 

OO 

^S"^ = E e'(''-''--(") with a(„),, 5,, - for all 



n. 



n=l 

oo 



^) = ^ 6(„), ,.(k.x-^<„,(k)*)^ 

oo 

j.("h) _ ^ ^^^^^ g»(k.x-<i(„)(k)t), (22) 



Above, 7r|^'' and z^^^^* have been expanded in the derivative expansion, which is an expansion in the scale of variation of 
hydrodynamic variables over the mean free path. We also require 5ui and 6T to be small uniformly for the linearized 
approximation to be valid. Furthermore, 77 is the shear viscosity and 2? is the charge diffusion constant. On the other 
hand a(„)jj, &(„)i and C(„) parametrize the dissipative non-hydrodynamic modes of the energy-momentum tensor and 
conserved current. The n here represents the various non-hydrodynamic branches of quasinormal mode perturbations 
which dissipate because their dispersion relations a;(„)(k), (D(„)(k) and W(„)(k) have negative imaginary parts. We 
require a(„)y /p, 6(„)i/p and C(„)/p to be small for the linearized approximation to be valid. 

We note the separation of iTij and Vi into hydrodynamic and non-hydrodynamic parts can also be done at the 

non-linear level. This is so because even at the non-linear level the hydrodynamic parts tt^^'' and vf^^ are solutions by 
themselves - from the perspective of kinetic theories this follows from existence of normal solutions as discussed before 
and from the point of view of gravity they give regular metrics via fluid/gravity correspondence. For any tt^ and 

the non-hydrodynamic parts t^^^^ and are just whatever remains after subtracting out the purely hydrodynamic 

parts TT^j^'' and vf^^ constructed algebraically from the profile of the hydrodynamic variables in the full solution of the 
energy-momentum tensor and conserved currents. 

In order to obtain the hydrodynamic modes at the linearized level, we simply put all a,(n)ij to zero in ( |2^ ) 

and impose the conservation of energy, momentum and charge : 

= 0, d^'j^ = 0. (23) 

We then obtain three modes, the sound mode, the shear mode and the charge diffusion mode. In the sound mode, 

5u(k) is parallel to k, 

^ \ d J e{T, p) + p{T, n) 



Se{k) =d5p{k)= ±Vd\6u{k)\ [e{T,p)+p{T,ti)j + 

Mk) ^ Sejk) 

P e{T, p) + p{T, n) 



(24) 



Above (...) refers to higher derivative corrections in powers of k. Using thermodynamic relations locally, one can 
obtain dT{k) and 5p{k) from Se{k) and Sp{k). 
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In the shear mode, 

i5u(k) is orthogonal to k, 
UJ = —i 



In the charge-diffusion mode 



Se{k) = 5p{k) = 6p{k) = 0. (25) 



Se{k) = 0, Sp{k) = 0, (5ui(k) = 0, 

cj = -iX>(T,^) I k |2 . (26) 



The quasinormal modes of the metric and gauge fields contains these hydrodynamic modes as the only branches 
in which uj and k can go simultaneously to zero. We can also obtain the transport coefficients by using the incoming 
boundary condition at the horizon. We will be interested in the shear mode in particular. The shear-viscosity is given 
by [|: 

s{T,fj) e{T,^l)+p{T,^i) An- ^ ^ 

Above, s is the entropy density and we have used the thermodynamic identity s = (e + p)/T. 

In order to obtain the sirnplest non-hydrodynamic modes we need to set the perturbations of the hydrodynamic 
variables Sui, ST and Sfi in (^2|) to zero. Also we look for spatially homogeneous perturbations so that the momentum 
k is zero. Nevertheless, unlike the case of hydrodynamic modes, the frequency do not vanish when k goes to 
zero. In such a configuration, for arbitrary a(„)y , it is easy to see that energy and momentum is conserved because 
d'^t^^ vanishes identically. When the chemical potential is set to zero, the quasi-normal modes in five dimensional 
gravity in AdS^ give |^ : 

a;(„) (k = 0) ^ ttT \± 1.2139 - 0.7775 i ± 2n{l T i)] , for large n. (28) 

Clearly, the conservation equations are not enough to reproduce all the quasi-normal modes. We need extra 
phenomenological equations. Such phenomenological equations can be derived from kinetic theories like Boltzmann 
equation at weak coupling or gravity at strong coupling. However, we can also write them on general phenomenological 
grounds. At present, these will not be important for us, we merely mention these have been found in the most general 
form in §, §. 

We will be interested in the spectral function in this class of non-equilibrium states, whose dynamics is determined by 
the non-equilibrium fluctuations of energy-momentum tensor and conserved currents only. If we want to obtain these 
spectral functions holographically, we need the explicit metric and gauge field corresponding to the non-equilibrium 
state. It will be important for us to write the metric and gauge field fluctuation about the equilibrium black-brane 

background explicitly in terms of Sui, 6T, J/i, t'o and t^l"'^'. As we will show in the next section, the spectral 

function in the dual states will depend explicitly just on these non-equilibrium variables. 

Later in section V, we will discuss what happens when we take into account non-linearities in the dynamics of 6ui, 
5T, TT^f^ etc. 



D. Explicit examples of backgrounds 

We will be interested in strongly coupled conformal field theories in three space-time dimensions in the large N 
limit. Therefore, as discussed earlier, we will be concerned with solutions of Einstein-Maxwell equations which are 
asymptotically AdS^ and are quasi-normal mode fluctuations about a Reissner-Nordstorm black brane with both mass 
and charge. 

As discussed earlier, on the gravity side we will need the Einstein-Maxwell action : 

^-2^/'^'K^+^-i^^^-^'")- ^'^^ 

Above I sets the scale of asymptotic (negative) curvature via a (negative) cosmological constant. This is required so 
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that the asymptotic isometry of the spacetime is the same as the conformal group in 3 dimensions. We will use k to 
denote the effective Newton's constant in four-dimensional gravity in lieu of Planck length Ip. 
The metric of the Reissner-Nordstorm black brane in AdSi is : 



where / is the so-called blackening function given by 



/(.) = i-^3;|+i).-v3;|a (31) 

^ 'o ' 'o 



In case of the gauge field, it is convenient to use the gauge Ar = Q. The only non-zero component of the gauge field 
is At and is given by : 



At = 



2v^ 



0-?)^ (32, 



The boundary of AdS^ in these coordinates is at r = and the outer horizon is at r = P /tq. The total mass M 
and charge Q of the black hole are given by : 

Q^VSrl M = rl + 3^. (33) 

ro 

Using the standard holographic dictionary we can relate the two parameters and tq of the geometry and the 
Newton's constant k in to the energy density e, charge density p and entropy density s as below : 



,3 , ^4 



e = 2p= -2^(3-4 +1), p^—{-), s=-^. (34) 

The thermodynamic relation 



de = Tds + pdp (35) 
gives the temperature and chemical potential as below : 



3ro / /'r^\\ 2^/3: 



6r 



2 



The first example of a non-equilibrium background we will describe is that with a hydrodynamic shear-mode turned 
on. The velocity perturbation will be denoted as (5u(/c(h)) with being the three- momentum of the fluctuation. We 
recall that k(h) ■ <5u(k(ij)) = 0, as the shear wave perturbation is transverse. 

It is a well-defined problem to find a given metric and gauge field perturbation in the bulk corresponding to a 
definite energy-momentum tensor and conserved current fluctuation about the equilibrium at the boundary, when 
the Dirichlet boudary condition is imposed for the bulk perturbations at the boundary. The latter is needed so that 
the dual field theory lives in flat space and is influenced by an externally fixed chemical potential. Regularity at the 
horizon fixes the transport coefficients appearing in the energy-momentum tensor and conserved currents. 

This procedure can be readily implemented in Fefferman- Graham coordinates |28[ . A similar procedure can be 
implemented in Schwarzchild-like coordinates as well because the Schwarzchild radial coordinate and the Fefferman- 
Graham radial coordinate are only functions of each other when the temperature remains unperturbed. Then it 
follows [in that : 

5gij will be proportional to (fc(ij)j (5uj(k(h)) -I- (5ui(k(ij)))e*^'^<''''''~"('')*-', and, 
5gM will be proportional to (5ui(k(h))e'(''('')'''""(")*'. 



It can be also shown that in the radial gauge, = 0, the fluctuation in the gauge field is also proportional to the 
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fluctuation in the conserved current, i.e. proportional to : 

(5M,(k(h))e*(''('"-^-'^<'"*). 

The exphcit metric is given by : 



+ ( " <57.,(k(h))e'(''<'"-"-"<'^)*' h (^) dx^dx^) + 0{e'), (37) 



where, 



h(s) = 3 / d§ 3 , (38) 

io (l + S + S2_3l|,54)(i_5)' 

and 

^(h) = -^^ + 0(.^), '? = ^=4-- (39) 
In the radial gauge Ar — 0, the gauge field takes the form 

= ^ (l - + Oie% = (l - 9^) ^..(V))e^(''<-x-.o.*) + Oie^). (40) 

Above e denotes the parameter of derivative expansion in hydrodynamics. 

It is to be noted that we have written the full metric and gauge field in a global frame co-moving with the equilibrium 
part of the energy-momentum tensor and conserved currents, i.e. in the laboratory frame. We can readily make the 
metric and gauge field Lorentz-covariant by boosting such that the unperturbed velocity field is a four- velocity vector 



Q. However, this will be unnecessary for the purposes of this paper as we will be interested in the results in the 
laboratory frame. 

Also one can readily realize that the metric is singular at the outer horizon r = P/tq. This is however only 
an artifact of the coordinate system. We can systematically change coordinates order by order in the derivative 



expansion so that the metric and gauge fields are manifestly regular at the horizon |28 . In our coordinates, the radius 
of convergence of the derivative expansion is of the order of the effective mean- free path or the inverse of the effective 
temperature at a given radius given by Teff{r) = T j \J /(rro/P). Therefore, we have a finite radius of convergence 
of the derivative expansion a finite distance away from the horizon. Furthermore, we will be interested in calculating 
boundary correlators which are independent of the choice of bulk coordinate system. 

The metric ( |37| ) and gauge field (|o|) in manifestly regular coordinates are given in appendix A. 

The second example which we will be concerned with will be a homogeneous non-hydrodynamic perturbation of 
the energy-momentum tensor, i.e. with one o,(n)ij in ( p2[ ) turned on. The momentum of this perturbation is zero on 
account of homogeneity, but its frequency is non-zero and complex like in ( |2^ ) . The metric can be obtained following 
[0 in the Fefferman-Graham coordinate and re-expressed in the Schwarzchild coordinate used here by simply changing 
the radial coordinate. Again, as the temperature remains unperturbed, up to linear order the change of coordinate 
involves transformation of one variable. It can be shown that the metric perturbation is proportional to 

Explicitly the perturbed metric is : 

P dr"^ P ! „ /rrn\ , , , o . n\ 21 



ds-" = 



^( - / (9^) dt' + dx' + dy^) _(a(„)„e--(")* h {^,u^in)) dx^dx^) + 0{6^), (41) 



with S being the parameter of non-hydrodynamic amplitude expansion. Furthermore, /i(s,cj(„)) follows the equation 
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of motion : 



d'h{s,U^^n)) (2 +(1 + 3^)^3- 6^.4) ^/,(,,^(^^^ 



|/i(s,L^(„)) = 0. (42) 



ds"^ sf{s) ds Tg I /^(s) 

We will also require that : 

/i(s,a;(„)) = / + 0(s*) ass^O. (43) 

This is the asymptotic boundary condition and determines h uniquely as it puts the coefficient of the non-normalizable 
to zero and the coefficient of the normalizable mode to be unity so that the boundary energy-momentum tensor 
fluctuation is as given by (p2|). Though the equation for h cannot be analytically solved, the solution can be readily 
expanded in a power series in cij(„) . 

Furthermore, the gauge field remains unperturbed from the black brane profile. 

The metric above is also not manifestly regular at the horizon, but once again it is just an artifact of the choice of 
coordinates. One can again translate the metric systematically to Eddington-Finkelstein coordinates to see manifest 
regularity |7|. The regularity is manifest only when we sum over all orders in a;(„). This is to be expected because, 
although the amplitude of the non-hydrodynamic perturbation a{n)ij is small, it's rate of change in time is not small 
(unlike the hydrodynamic modes) since i^(n) is of the same order as the temperature. 

Though we will not discuss the details here, we can construct the explicit metrics in the case of both hydrodynamic 
and non-hydrodynamic perturbations even at the non-linear level 0|. The metric is regular at each order in 
the derivative expansion for hydrodynamic perturbations and for each order in the amplitude expansion for non- 
hydrodynamic perturbations, provided all time-derivatives (or covariantly speaking convective derivatives) are summed 
over at each order in the latter case [fzl . 



III. THE HOLOGRAPHIC PRESCRIPTION FOR THE NON-EQUILIBRIUM SPECTRAL FUNCTION 

As discussed in the Introduction, the spectral function is given by the imaginary part of the retarded propagator 
which can be obtained from causal response of an operator to it's source. A convenient way to obtain the spectral 
function is to calculate the retarded propagator using linear response theory first and then isolate its imaginary part. 

In this section, we will consider single trace scalar and fermionic operators in field theory whose back-reaction to 
the metric is suppressed by 0{1/N'^). As we have argued in section II. B, the possibly interesting scalar operators in 
the strong coupling and large N limit are order parameters of symmetry breaking. If we are in a range of temperature 
and chemical potentials, where such symmetry breaking does not occur, the profile of the scalar fields dual to these 
operators vanishes in the background classically. Therefore, the backreaction is indeed 0{1/N'^) suppressed. This 
observation may be applied to study pion correlations in the quark-gluon plasma at RHIC. 

In popular holographic models of strongly correlated systems, the electron is thought to couple to a composite 
operator made out of strongly interacting fractionalized degrees of freedom (for a clear exposition please see p3| ). 
The holographic dual is thought to capture the dynamics of the fractionalized degrees of freedom. The strongly 
interacting fractionalized degrees of freedom are 0{N^), but the coupling of the electron to the composite operator of 
the strongly coupled theory is 0(1). The spectral function obtained from photo-electron spectroscopy (ARPES) will 
receive corrections from the spectral function of the composite fermionic operator of the strongly coupled sector. As 
the coupling of the electron to this operator is 0(1), we can ignore the backreaction of the fermionic field dual to this 
operator on the geometry representing the dual state, at the leading order. If this picture is qualitatively viable, our 
set-up will be relevant for describing non-equilibrium features of non-Fermi liquids described by such models. 

Holographically, causal response implies the incoming boundary condition at the horizon. The event horizon 
separates space-time into two causal parts, one that is inside and ends at a singularity, and the other that is outside 
and stretches all the way to the boundary. No light ray can come out of the inside region to the outside region, 
though light rays can propagate from the outside to the inside. Therefore, the perturbations which respect the causal 
structure of the space-time are those which are purely incoming at the horizon, having no component which propagates 
from the inside to the outside. 

The event horizon is not only a feature of the eternal static black hole, but also of the perturbed black hole (for 
instance, the black hole with the quasi- normal mode fluctuations of the metric and gauge fields). The event horizons 
of these non-equilibrium geometries are also perturbed from their equilibrium location and their positions can be 
calculated in a perturbative expansion p^ ]. Equilibration in this context means that the event horizon will have 
uniform surface gravity (the gravitational analogue of temperature) everywhere and it happens only far in the future. 
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Though the incoming boundary condition is insufficient for a well defined perturbation theory in non-equilibrium 
geometries as noted in the Introduction, we expect regularity at the future horizon to be a sufficient condition. It 
turns out that it is sufficient to impose the regularity condition only far in the future, that is in the asymptotic static 
black brane geometry. This has been observed before in |^ in another context - while constructing time-dependent 
non-linear solutions of gravity with regular future horizons perturbatively. In such solutions it indeed suffices to 
impose regularity of the perturbations at the final equilibrium location of the horizon. In fact, the incoming boundary 
condition is itself tied up to regularity p5[ ^ . In this section we will find a precise non-equilibrium generalization of 
the incoming boundary condition for bosonic and fermionic field configuratons in non-equilibrium geometries. 

For purposes of illustration, let us consider the non-equilibrium state which is the simplest to analyze from the 
gravity point of view - it is the AdS black brane with a linearized hydrodynamic shear mode perturbation of spatial 
momentum k(jj) . The advantage of this geometry is that it can be shown that the event horizon do not fluctuate up 
to first order in the derivative expansion (i.e. up to first order in k(h)/T) essentially because the temperature field 
does not fluctuate as discussed in section II. We will flrst demonstrate how we can develop a prescription for obtaining 
the holographic spectral function in such a non-equilibrium state. Our aim will be to obtain the correction to the 
equilibrium spectral function up to first order in derivative expansion, i.e. up to first order in k(ij)/r. 

The explicit metric and gauge field of the black brane with the hydrodynamic shear mode perturbation is given in 
(^) and ( ^o| ) respctively up to first order in the derivative expansion. We will work explicitly with four space-time 
dimensions in gravity, as we will be interested primarily in a three space-time dimensional dual strongly coupled field 
theory. This is because we are interested in applications to strongly correlated electron systems at finite density 
living in two spatial dimensions. As argued in section II. B, our analysis may apply to the strange metallic phase in a 
qualitative manner. 

An elegant way to solve the equations of motion of scalar and fermionic fields is by using the Fourier transform in 
all the field-theory (i.e. boundary) coordinates. Obviously, in order to express the equations of motion of the fields 
in Fourier space, it is necessary to do the Fourier transform of the background perturbation first, i.e. we need to do 
the Fourier transform of the velocity field fluctuation Sui . The dispersion relation for this fluctuation is as given by 
eqs. ( psi) and (|2^). We see that the frequency given by the dispersion relation is strictly (negative) imaginary, while 
the frequency related to Fourier transform is strictly real. Furthermore, the negative imaginary frequency given by 
the dispersion relation makes Sui decay in the future but grow in the past as a function of time. A Fourier transform 
of such a function needs to be defined with care. In order to distinguish from the frequency and momenta associated 
with the scalar/fermionic field, we will denote the frequency and momenta of Sui as w^h) and k(i,-) respectively. The 
correct Fourier transform which reproduces the hydrodynamic dispersion relation is : 



To check the above, one can try to reproduce the time dependence by doing the inverse Fourier transform. This needs 
to be done with a specific contour prescription for integration over a;(h) a-s shown in Fig.0.This contour is the usual 
contour associated with the retarded propagator in field theory - it runs from — oo to oo infinitesimally below the real 
axis and then closes itself through the circle at infinity. This contour picks up contribution only from the negative 
imaginary pole reproducing the correct time dependence of Sui at given k(h) • 

It will be easier to solve the scalar/fermionic field equations after doing the Fourier transform of Sui, however we 
need to finally integrate over tLi(h) with the above contour prescription in order to obtain the observed behavior in 
real time. 

For demonstrative purposes, we will analyze the scalar field equations first and then the fermionic field equations. 
Finally, we will see how we can apply our prescription for the non-equilibrium retarded Green's function when the 
background contains other quasinormal modes of the metric and gauge field. 



We will be interested in the non-equilibrium holographic spectral function for a scalar operator first. This requires 
us to solve the equation of motion of the dual scalar field in the non-equilibrium background; in particular we need 
to understand how the equilibrium part determines the non-equilibrium part completely. Without this, as we have 



See also M for an explicit proof in a non-hydrodynamic context. 




(44) 



A. Scalar field equation and the non-equilibrium spectral function 
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FIG. 1: Contour for integration over i^(h), with pole at negative imaginary axis 



mentioned before in the Introduction, the spectral function cannot be determined. 

We will need to specify the equilibrium part of the solution first. We can assume, without loss of generality, that the 
equilibrium solution is in a specific {uj, k) mode and obtain the non-equilibrium correction for each such mode. Using 
the fact that our field equation is linear, we can then linearly superimpose the solutions with the non-equilibrium 
correction for each equilibrium mode to obtain the most general solution. 

The background in which the scalar field propagates is the AdS/i Reissner-Nordstorm black hole with the hydro- 
dynamic shear-mode perturbation. This hydrodynamic mode is given by the velocity perturbation 5ui in a specific 
momentum k(h) but its dependence on w^^) is given by (Q). We have to consider the background first in a definite 
perturbation and then integrate over uj^y^ finally with the contour prescription discussed before. The scalar field 
while propagating in the background will pick up a (oj -f W(ii),k + k(h)) mode. The profile of the scalar field, will 
therefore be of the following form : 



$(x,t,r) = $("'(a;,k,r)e-'("*-''-^) -h$(^Ha;,k,tj(h),k(h),r)e-*(('^+"(''))*-(*'+''<>'))-^^ 



(45) 



The equilibium part of the solution is $('^)(cj, k, r) and the non-equilibrium part is $'^^''(w,k, cij(h)ik(h),r). The non- 
equilibrium part does not depend on the combination ui + a;(h) and k -I- k(h) as the space-time translational invariances 
of the equilibrium background are broken explicitly by the hydrodynamic quasinormal modes. 

If the scalar field $ is minimally coupled to gravity, and its mass and charge are m and q respectively, the equation 
of motion of the equilibrium part is simply 



Utf'^,5{u' - u)5\\^ - k)$(°) (a., k, r) = 0, 



(46) 



where D^k^ is the (gauge-invariant) Laplacian in the AdS Reissner-Nordstorm background metric (30) and gauge 
field (32) as given by : 



□ 



ARN 



dr 



2 i2 

r I 



(a; + q^i (l 



rra\\2 



(47) 



Again, / is the blackening function of the AdS Reissner-Nordstorm black brane which vanishes at the horizon located 
at r = P /tq. 

With the metric and gauge field in presence of hydrodynamic shear perturbation given by (^^ and ( ^0[ ) respectively, 
the equation of motion for the non-equilibrium part up to first order in the hydrodynamic momenta '^(h) is : 
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^t^^' - ^ - ^(h) (k' - k - k(h) )<I>(i) {io, , k, k(i,) , r) = ^(w, 6.(1,) , k, k(h) , r)<i>(") (l., k, r) , (48) 

with 

V = ^1+^2, 

1^1 = 



(^(l - ^(9^)) + '^'^ " 9^) )Mw(h),k(h)) • k, 

^;:;^^(72~j^^%^(h)i'^"i('^(h)>k(h))- (49) 

Above, /i gives the hydrodynamic correction to the background metric which is proportional to k^Y^-^^5uj + («<->■ j) as 



in (|38|). 

The behavior of the general solution of $^°'(a;, k, r) near the horizon is well-known. It can be split into an incoming 
and outgoing wave as below : 



$(°)(w,k,r)«A'"(a;,k)|^l-^j + A°"*(^, k) |^1 - ^ j near r = ^. (50) 

In order to select the incoming wave, we should put 

A°"*(w,k)=0. (51) 
We can also normalize the overall solution by choosing 

^™(c.,k)=C, (52) 

with C being a numerical constant. This overall normalization will play no role in the Green's functions. 
The behavior of the general non-equilibrium part of the solution near the horizon is : 



;2 



"(h) 




$W(u;,u;(h),k,k(h),r) « A™(u;, k, k(h)) [ 1 - ^ j + A°«'(c., k, k(h)) (^1 - ^ j 

\ w(5u(w,k(h)) -k / rro\ '""^ 

1 -5- nearr = — . (53) 

(2a; + a;(h) )t.(h) \ ' / '"0 

The first two terms on the RHS above are the homogeneous incoming and outgoing solutions for frequency mode 
Lo -\- ■ The third term is the particular solution which is determined completely by the equilibrium solution. The 
above behavior at the horizon is exact up to first order in k(i,) . In fact the full general solution which reproduces the 
above can be given elegantly in an integral representation as in appendix B. 

Obviously, we need to impose the incoming boundary condition again. Therefore, 

A°«*(c.,c.(h),k,k(i,)) = 0. (54) 

We will now show that in order to impose regularity at the horizon, we also need to dispose of the ingoing non- 
equilibrium homogeneous solution at the horizon. We recall that finally we need to integrate over a;(h) . 

In order to be consistent with the derivative expansion, A*"(w, ^(h-), k, kjh)) must take the form as follows. It is 
proportional to components of i5u at the linear order as it should vanish in absence of the background perturbation. 
It's dependence on a;(h) and k(h) can be expanded systematically in terms of rotationally invariant scalars like Su ■ k, 
kikjk(yji6uj, oj (\y)kikjk(yiiSuj , etc. Up to first order in the derivative expansions only the first two scalars will apear. 
The coefficients of these scalars should be functions of lo and k only, as the depenedence on a;(h) and k(h) can be 
absorbed in coefficients of the scalars appearing at higher orders in the derivative expansion. Thus, up to first order 
in derivative expansion, we should have : 
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^"(w,W(h),k,k(h)) =ylf(cj,k) (5u(a;(h),k(h)) -11 + A™(a;,k) fc,fcjfc(h)i(5uj(w(h), k(h)). (55) 

We recall for the hydrodynamic shear mode Su ■ k(ii) = 0, so there is no more possible terms up to first order in k(ii). 
When we integrate over a;(h) , the Fourier transform of Su as given by ( ^ ) will give a pole contribution. Taking this 
into account the behavior of the ingoing non-equilibrium mode at the horizon will be : 



(56) 

Therefore, we find the ingoing homogeneous non-equilibrium mode diverges at the horizon as k^j^^/ (IGtt^T^) is strictly 
positive. This divergence is not an artifact of the coordinate system because we are studying the behavior of a scalar 
field. The only way this divergence can be removed is by putting 

A'"(^,^(h),k,k(h)) -0, i.e. Ar(c.,k) = A™(c.,k) =0. (57) 

The particular solution at the horizon as defined as the third term in (^3|) produces no divergence after we do the 
integral over a;(h) . It is regular at and outside the horizon. 

Summing up, the full solution with the non-equilibrium correction is the following : 



H.,t,r)^C 1-^-^] ^^-.--)+. 2 WM.,k,,).k/ ^ 

g-»((.^+a;(^))t-(k+k(^)).x) j ^ ^g^j. ^ = (58) 

J ro 

The above behavior when specified near the horizon uniquely fixes the full non-equilibrium solution aside for an overall 
normalization C. 

We can numerically extrapolate the full solution all the way to the boundary r = 0. As the background is 
asymptotically AdS, we should have the following behavior : 

$(x,i,r) « J(x,t)r^"^ + 0(x,t)r^ near r = 0. (59) 

By the holographic dictionary, J is indeed the source and O is the expectation value of the dual operator in the dual 
non-equilibrium state ^. Also, A is the scaling dimension of the dual operator given by the mass of the scalar field as 
below : 



A = ^ + ^^+m^P. (60) 

The positivity of the Hamiltonian requires m^/^ > —9/4 

Furthermore, near r = 0, the equilibrium and non-equilibrium parts of the solution individually have the same 
behavior, so 

$(")(tj,k,r) w J(°)(cj,k)r3-'^ -|-O(0)(tj,k)r^, 
$(i)(L.,L.(h),k,k(h),r) « j(i)(L.,c^(h),k,k(h))r3-^+0(iHc.,c.(h),k,k(h))r'^. (61) 



When -9/4 < ■m?l'^ < -5/4, we can do an alternate quantization where J can be interpreted as the expectation value and O as the 
source |36[ . This requires the scaling dimension of the operator to be A = 3/2 — 1^9/4-1- va^l'^. The partition functions of the two 
theories are related by a Legendre transform. 
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Therefore, 



-/ — OO 

0(x,t) = 0(")(Lj,k)e-*('^*-''-") + / dw(h)0(iHcj,cj(h),k,k(h))e-*«'^+"(>'))*-(''+''(>'))-'^). (62) 



The unique solution of ^^^^ with our prescribed behavior near the horizon (^8|) gives us the precise non-equihbrium 
contributions to both the operator and the source in the fohowing form : 

0^^)(a;,a;(h),k,k(h)) = o'x^(u,]ij 5u(a;(h), k(h)) -k + o'-g^ (^uj.kj fcjfcjfc(h),,5uj(a;(h), k(h)), 
J^^)(a;,a;(h),k,k(h)) = J^^ (^w, k^ (5u(w(h), k(h)) • k + J^^^ (^w, k^ fci%fc(h)i(5wj(a;(h), k(h)). (63) 

The exphcit forms of o'"^\ Og\ jjj^^* and jjj^-* can be obtained as in appendix B. The integration over W(h) then will 
be given by the contribution from the pole in 6u. 
The non-equilibrium retarded correlator is 



G7?(xi,ti,X2,i2) 



where 



0(xi,^i) 



-iw(ti-t2) ik-(xi-X2) 



0(°)(w,k) + 0(i)(^>k,k(h))e*<")-^ie-^*i 

P 

J(")(a;,k) + J(i)(tj,k,k(h))e*(")-''2g- 4^-42 



g-iLj(t i-t2)gik-(xi-X2) 



O(o)(cj,k) 



J(o)(w,k) 



'0(i)(^,k,k(h)) 
0(0) (w,k) 



(h)J „»k„vxi 



(ii) , 



j(^)(^,k,k(,)) , 

J(0)(a;,k) 



(64) 



0(i)(c^,k,k(h)) = 0i'^(c^,k)5u(k(i,)) -k + 0i'^(c<.,k) fc,fc,fc(h).(5u,(k(h)), 

j(^)(a;,k,k(h)) = jj^^^ (^w, k^ 5u(k(h)) • k + (^w, k^ fcjfcjfc(h)j5wj (k(h)). (65) 

The difference of the above from (^3|) is that in 5u which has no dependence in a;(h) ■ The latter has been integrated 
over. This integration produces the contribution from the diffusion pole and the residue has been obtained from (^). 

Clearly, the choice of overall normalization of the solution given by C in (^8|) does not matter as mentioned before. 
It cancels between the numerator and denominator in the retarded correlator. To readily compare with experimental 
data, we have to do the Wigner transform of the retarded correlator, as discussed before. We find 



G'fl(i:j,k,x,t) = j duo j (fka 



J(°H^o,ko) 



OW(L.o,ko) 1 / C>(i)(a.o,ko,k(h)) ,/ k(h) 
" j(o)(a;o,ko) 27rzl^ O(o)(wo,ko) V ° 2 



SttT 



jW(^oko,k(,)) . 1^ 
J(o)(c<.o,ko) V 2 



|g»k(h)-Xg-^t 



(66) 



The first term above is just the equilibrium retarded propagator. The second and third terms are the non-equilibrium 



At equilibrium, this prescription has been proposed in p^ . As noted in the Introduction, we can apply this prescription also at 
non-equilibrium using the validity of linear response theory. 
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contributions. The non-equihbrium contributions have an exphcit space-time dependence which is co-moving with the 
velocity perturbation in the background. 

The spectral function can be obtained from the imaginary part of the retarded propagator by using A{uj, k, x, t) — 
-2ImGfl,(w,k,x,i). 



B. Fermionic field equations and tiie non-equilibrium spectral function 

We will now extend the prescription to obtain the non-equilibrium fermionic spectral function. We begin by 
constructing the equation of motion for a Dirac spinor explicitly in the same non-equilibrium background, which is 
^^5*4 Reissner-Nordstorm black hole with a hydrodynamic shear-mode perturbation. 

We recall that the Dirac equation for a Dirac spinor of mass m and charge q in curved space is : 

e^'r^ (dM + lujff[rB,rc] + iqAm] + m] * = 0, (67) 



where M are the space-time indices, and A, B and C are the tangent space indices collectively. We will denote tangent 
space indices with underlines as in (r , i , a;, y) or more compactly as (r, /i) to distinguish from the space-time indices 
which will not be underlined as in (r, t, x, y) or (r, ^). 

In order to work with the holographic dictionary, it is convenient to choose the following representation for Gamma 
matrices [ p3[ : 

where 7^s are the 2-1-1 dimensional Gamma mtrices in a chosen representation. We will choose the latter in the 
following representation : 

It is also useful to decompose the 3-1-1 space-time dimensional Dirac spinor as eigenvectors of T± defined as : 

r± = i(i±r^), (70) 

so that 

* = + *±=r±*. (71) 

The advantage of this decomposition is that both ^1*+ and transform as 3 space-time dimensional Dirac spinors 
when the Gamma matrices are in the representation above. 

It might be puzzling as to how a Dirac spinor in the bulk maps to two Dirac spinors in the boundary, but we note 
unlike the scalar field equation, the Dirac equation is first order. Therefore, as in the case of the scalar field we have 
two independent boundary data, corresponding to 5*+ and each. Eventually, we will see how these two boundary 
data maps to source and expectation value of the dual operator, and further how they get related to each other by 
regularity in the bulk giving us the dual fermionic retarded propagator. 

Just as in the case of the scalar field, the space-time profile of the Dirac spinor also has an equilibrium and non- 
equilibrium part. We can first assume that the equilibrium part is in a specific {u, k) mode and determine the 
non-equilibrium correction to this. Later, we can obtain the most general solution by superimposing the full solutions 
corresponding to various equilibrium modes. The space-time profile of the Dirac spinor thus takes the following form 



*(x,t,r) =*(*')(w,k,r)e"*("*-''-") -|-*(i)(w,k,W(h),k(h),r)e-'«"+'"("))*-(''+*'(^))-"\ (72) 

where '^^^^ is the equilibrium part, Vf^^^ is the non-equilibrium part, and (w(h),k(h)) correspond to the frequency and 
momenta of the velocity field perturbation in the background. From now on, we will denote (w,k) collectively as fc, 
and (aj(h)i kj-ii)) collectively as fc^jj). 

The equations of motion for can be written as two coupled first order PDEs for It will be convenient for us 
to decouple these PDEs and write a second order PDE for 4'+. It will turn out that will be then algebraically 
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determined by For the equilibrium AdS4 Reissner-Nordstorm black brane background, this has been done in 
[Q. Following this, we write the equations of motion for as below : 



-A" 



vff (fc,r), 



(73) 



where 



Qik,r) 



A+A- 



ro 



A^ 



k 



(74) 



and 



A" 



tk 



1 

27 



± 



=m. 



(75) 



with ' denoting differentiation w.r.t. rrp/P, representing the equilibrium configuration of the gauge field and 7^ 
is tktk- 

In order to obtain the equations of motion for ^I'^'' we need to obtain the non-equilibrium first order corrections 
to the vielbeins and spin connections in the derivative expansion. These are given in details in appendix C with the 
metric being ( ^7| ) corresponding to the black brane perturbed by the hydrodynamic shear mode. 

In order to simplify calculations, we will choose (without losing any generality) the momentum of the velocity field 
perturbation k(h) in the background to be in the x direction; therefore the velocity perturbation 5\i being transverse 
should then be in the y direction. Later, we can make the results manifestly rotationally covariant by rotating, and 
also Lorentz covariant by boosting to an arbitrary frame. The momentum of the equilibrium part of ^ of course can 
have arbitrary components in both x and y directions if we have to retain full generality. 

The equations of motion of are as follows : 



^ + Pik',r)^ + Q{k\r)]s^k' -k- fcy,)) *W(fc, r) 




5+(fc, fc(h),r) 



5a(fc'-fc-fc(j^))vI/W(fc,fc(h),r) = f|;+^+U3(fc'-fc-A:(h)) vI'V^(fc,fc(h),r) 



S+{k,k(^i,^,r), 



(76) 



where 



5+(fc,fc(h),r) - -A'+(fc(h),r)vI/f (fc,r)-3;(fc(h),r)*^'^(fc,r) 
5_(fc,fc(h),r) = -A'_(fc(h),r)vI/L°^(fc,r)-3^(fc(h),r)vI.f (fc,r) 



(77) 
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with 



A'±(fc(h),r) = Ti(f(A;(h),r)7^7^--^(fc(h),r)7^^), 
yikii^),r) = ^(s(A:(h),r-)7^-C(A:(h),r)7^7^^) +e(fc(h),r)7^ + H(fc(h),r)7^. (78) 
B, C, £ , T ^ QjH are given in terms of the inverse vielbeins and the spin connections as 

I3ik(h),r) = — ^ (2iqelAy + 2iqelAt+elujf- + 2iujel-2ik^ef'^^^ 
'S^(fc(h),?-) = — , [~^t^T ~e«y^^ + er^r--eyu;r-et^y ) 



r 



rrg 



^ fx y— y xr X xr i v y—\ ^ ^ 



/(^) 



?^(fc(h),r) = J (-»fc,e|)^'^ (79) 



Here (• • • )^^-' means that we are extracting only those parts of the full expression which is first order (i.e. linear) in 
k(h). Once again we mention that the exact expressions of the inverse vielbeins (or einbeins) and spin connections 
appearing above are given in appendix C exactly up to first order in k(-h) . 

The most important observation regarding the equation of motion for ^'^^^ is that just as in the case of 
as evident from (^6|), can be determined first by solving a second order ODE and can be determined 

algebraically in terms of the solution for ^_|_. Therefore to uniquely specify ^'^^ it is sufficient to uniquely specify 
^I'^''. Moreover, the differential operator on the LHS of the equation of motion (|7^) for is the same as that for 
^'^^ in ( [73I ) with k replaced by fc + fc(h)- Therefore, the homogeneous solutions of ^'"^^ will be the same as those of 
^''^^ with k replaced by fc + fc(h) . 

The general behavior of the equilibrium part of the solution '^^^^ at the horizon r — l'^ /tq is 

^^\.,Kr) « A!r(.,k)(l-^)"'^"^Ar(^,k)(l-!^)'^-^ (80) 
Both A™ and A™* are arbitrary linear combinations of 

0) (1 

The incoming wave boundary condition requires us to impose 

A™*(w,k)=0. (81) 
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Furthermore, the choice of k) will not matter in the final answer for the retarded propagator, so we will choose 

A-(c.,k)= (^^) (82) 

with /C being a constant. The behavior of near the horizon can be obtained via the second algebraic equation of 
( [73I ) as below : 

^^\.,Kr)^-,^-{l-'-;^y^-''[^). (83) 

Thus ^''l'^ is also incoming at the horizon and is 'i!^^'^ times a specfic function of the frequency and momenta. 

It is to be noted that the incoming wave solution of the fermion diverges at the horizon as well. That this divergence 
is not an artifact of choice of coordinates can be seen by computing the scalar 'I'vE' at the horizon. In fact, it is believed 
that the fermion backreaction at the horizon is strong enough to change the near horizon geometry of the black brane 
p8[ . As mentioned in the beginning of this section, we will assume here that the backreaction is suppressed by a 
factor of 0(1/7V2) 9 



We now turn our attention to the non-equilibrium part of the solution. From, the first equation in (76) we obtain 



that near the horizon, ^V^'' behaves as 



vI/W(c.,k,L.(h),k(h),r) « A™(c.,k,c.(h),k(h))(l-^) '+^r(w,k,L.(h),k(h))(l-^) 

+a(a;,k,tj(h),k(h))(l- ^) ' (^q^ 



a(.,k,.(,,k(,,) . 8(3.r + 4(7.T;..) ^^-^(^C^))^^^- («4) 

when we have chosen the incoming wave boundary condition and our normalization for "^^^K Thus we have again two 
arbitrary coefficients for the incoming and outgoing homogeneous solutions at 2 + 1 momenta k + j and then we 
have a particular solution completely determined by the source term. 

We now apply a similar logic as in the case of the scalar field. We put in (|4[) to be zero again to satisfy the 
incoming boundary condition. In order to be consistent with the derivative expansion, j4*" has to linear combinations 
of (Su(a;(ii), k(h)) • k and kikjk(y^i5uj{uj(\i^^ ''(h)) with coefficients which are functions of w and k only. The integration 
over a;(h) in presence of 5\i will give contribution from the diffusion pole which will cause a further singularity in the 
behavior of the fermionic field. This singularity involves an extra factor of 



So we put A™ to be zero too. There is however, a difference in the behavior of the particular solution near the horizon 
from the scalar evident from (IsJ) . It diverges at the horizon with an extra factor of 



The situation, therefore admittedly is confusing as both the incoming homogeneous solution and the particular solution 
are divergent by an extra power. Moreover, for sufficiently small hydrodynamic momenta k(i,-|, the divergence of the 
particular solution leads over that of the incoming homogeneous solution. 



At order 0(1/ N'^) we cannot ignore the backreaction even in the hnearized hmit. This is because the scalar and fermionic fields have 
non-trivial profiles even in the background due to Hawking radiation. Particularly, the Hawking radiated fermions forms a Fermi-sea in 
the near-horizon region of the AdS Reissner-Nordstorm black hole. 
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Nevertheless, we can argue as follows. When we take the backreaction into account, the part of the non-equilibrium 
solution completely determined by the source can be expected to be regular, as the source involving the regular 
equilibrium solution in the modified background will be regular in the next order in perturbation. This feature is 
observed in the case of fluid/gravity correspondence or for more general time-dependent solutions in gravity - if we 
make the solution regular up to 71— th order in perturbation theory, the source terms in the equations for n + 1— th 
order perturbations are also regular, and the divergences at the n + 1— th order can be removed by adjusting the 
homogeneous solutions only |2^ . 

In the present case, we will argue that the divergence of the incoming homogeneous piece coming from the integration 
over Ld(i^-f is there as long as the backreacted background has a horizon at the zeroth order. If indeed there is a 
horizon, we can define an incoming wave also through geometrical optics approximation. We can certainly construct 
an appropriate function of r which we denote as r«(r) such that the incoming radial null geodesic at the (modified) 
horizon is : 

V = t ~ r^{r). 

Clearly (r) has to increase indefinitely as r moves towards the horizon because of blue-shifting. The incoming wave 
at the horizon will always behave like : 

g-*("+"(ii))'' 

as the geometrical optics approximation is always good at the horizon due to the blue-shifting. Therefore, as long as 
the backreacted geometry still has a horizon, the integration over w^^-j will produce a divergent factor : 

(h) 

(r*(r)) i6„-^T-^ . 

Above we have used the result that the hydrodynamic dispersion relation up to the leading order remains the same 
in the presence of backreaction as "q/s is universally l/47r in Einstein's gravity minimally coupled to any form of 
matter Therefore, this divergence is not removable by backreaction as long as we do not get rid of the horizon 
completely. 

Getting rid of the horizon is generically impossible if we demand that the solution in gravity is well behaved, as 
that would expose the singularity unless the latter is also removed by the backreaction. The removal of singularity by 
back-reaction is impossible in Einstein's gravity minimally coupled to well-behaved matter. It is also hard to argue 
that solutions in gravity with naked singularities could be dual to states in thermal and chemical equilibrium in the 
dual theory. 

We conclude that the sensible thing to do is to proceed as in the case of the scalar field and put both A™ and 
to zero in the non-equilibrium part of the solution. This determines ^P^j^' completely and its behavior near the horizon 
is : 

(1)..,. V,., .^ - A^ K)(3z^T-^ + ^(,))-2(7r^r^+^^)) 



^'i'Va;,k,Lj.M,kfM,?-) = J^^-^ ^ ^ ^JuJfcfM)7^- 

+ ^ ' ' ^^>' ^^'' ' ^1 ttT 8 {SttT + ioj) {7ttT + iuj) ^'^>" ' 

1 — -j^^ ""^ " ^g'^ + sub- leading terms. (85) 

Once Vl'^'' is completely specified as above, we can determine ^'L^'' readily from the second equation in (^^ as it is 
algebraic. The behavior near the horizon is given by : 

*L^^(a;,k,W(h),k(h),r) = 



(1)^ ^ r77(27riT-2w + a;(h)) (igTT^T^-f lli7rTw-2w2-ha;(h)(2i7rT-w)) 

'■' '■'(h),k(h)''^) - 8 (SttT + luj) (TttT + iuj) (w + W(h)) 



5uy(fc(h))7-(^l — * + sub-leading terms. (86) 

We can integrate numerically from the horizon and find the full profile of (both equilibrium and non-equilibrium 
parts included) all the way up to the boundary. 

At the boundary, the behavior of is specified completely by the AdS^ asymptotic nature of the background. 
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When TO > 0, the behavior of 4'+ at the boundary is : 

^+{k, k^^^,r) « (j(°)(fc) + J^^Hk, fc(h)))r3-^ + (xW(fc) + M(i)(fc, fc(h)))r^+\ (87) 

with A being the scahng dimension of the dual operator and is related to the mass of the fermionic field by : 

3 

A = -+ml. (88) 

Clearly and A^^^) are determined by and J*^^-' and A^^^^ are determined by Similarly, the behavior 

of ^t- at the boundary for to > and to ^ 1/2Z is : 

vl/_(fc,fc(h),r) « (AA(0)(fc) +AAW(fc,A:(h)))r4-^ + (o^"' (fc) + (^^ fc(h))) r^. (89) 

When to = 1/2^, the leading powers of the homogeneous solutions above become the same. The behavior of at 
the boundary is then given by : 

(A:, fc(h) , r) « (aA(°) (k) + AA(i) (fc, A:(h))) r^ln r + (o^°^ (fc) + O^^) (A:, fc(h))) r^. (90) 

As is determined by '5+ algebraically, we get 

0(fc,fc(h)) = -^(2TO + l)X(fc,fc(h)), AA(fc,fc(h)) = ^2!^^J(fc,fc(h)), j-k = j^k^, e^k'^k^, (91) 

where O = 0(°) +0(i), etc. Thus we have just two independent boundary data corresponding to the fermionic source 
and expectation value of the fermionic operator dual to the field. The holographic dictionary indeed identifies J as 
the source and O as the expectation value of the operator when m > . Both these are fixed up to an overall 
normalization constant by the incoming boundary condition at the horizon and our regularity argument. 

Changing the sign of m is equivalent to interchanging 4'+ with |2^. Consequently J gets interchanged with O, 
and Ai gets interchanged with Af When to < 0, the scaling dimension of the dual operator is given by : 

A = ^-ml. (92) 

Once the solution in the bulk is determined, the source and the expectation value of the fermionic operator get 
related by a matrix V : 



J(^a;,k,W(h),k(h)j = I^l^a;, k, k(h) j 0(^a;, k, k(h) j . (93) 

Clearly D is independent of the choice of A"' for the equilibrium solution as we have claimed earlier. The retarded 
propagator is given by : 

Gfl(^w,k,a;(h),k(h)^ = i2?(^w, k, a;(h), k(ii)^7-. (94) 

Furthermore, as the non-equilibrium part of the solution is completely determined by the equilibrium part of the 
solution, we can compute the relations : 

0(^^(^w,k,cj(h),k(h)) = 7^A(a;,W(h),k,k(h)^ 0'-°\uj,k), J^^^ (^w, k, k(h)) = (t^, t^(h), k, k(h)) j'-°\uj,k). 

(95) 



When < \m\ < 1/21 we can also do an alternate quantization in which O is interpreted as the source and J as the expectation value. 
This requires the scaling dimension of the dual fermionic operator to be A = 3/2 — |m|Z. The partition functions of the two theories are 
related by a Legendre transform. 



28 



Above TZa and TZb are fully determined by our boundary conditions on ^J''^'^-' at the horizon. They take the form : 

7?.A('^,'^(h),k, k(h)^ = 7?.AA(w,k^(5u(w(ii),kh) • k + 7^AB (^w, k^ fcifcjfc(h)i(5-Uj (^(h), k(h)), 
TZb [uj, uj{h) , k, k(h) ) = T^ba (t^, k) (5u(tj(h) , kh) • k + 7?.bb (w, kj kikjkt^Yi)i5uj (w(h) , k(h)). 



(96) 

By going through the steps as in the case of the scalar field, we can easily see that the generalization of the form 
of the bosonic non-equilibrium retarded propagator (pq) to the fermionic case is : 



where 



(wo, ko)7-<5(w - a;o)5'(k - ko) 
- ^ (wo , ko )7^7^A (wo , ko , k(h) ) (k - ko - ^ 

-Ub (wo, ko, k(h))l?(°) (wo, ko)7-<5' (k - ko + ^) 

7^A(^^,k,k(ll)^ = 7^AA(w,k^(5u(kh) • k + 7^AB(t^,k^fc^fcifc(h)^'5MJ(k(h)), 
7^B(w,k,k(h)) = 7^BA(t^,k)^u(kh) • k + 7^BB(w,k)fcjfcjfc(h)i(5uj(k(h)). 




(97) 



(98) 

The first line in (^7|) denotes the equilibrium correlator and the lines below are the non-equilibrium contributions 
co-moving with the background velocity perturbation. The difference between ( p8| ) and ( p^ ) is that the integration 
over W(h) has kept only the residue of the diffusion pole in the Fourier transform of 6\i given by (^). 

Once again the spectral function can be obtained by computing the imaginary part of the retarded propagator 

above and using A(w,k, x, t) = — 2Im(Tr(7-Gfl,(w, k, x, t)) 



C. Generalization to backgrounds with other quasinormal modes 

The prescription we have presented so far is for the non-equilibrium retarded propagator in the hydrodynamic 
shear-wave background. We will now show that this prescription with it's underlying logic can be readily generalized 
to any background which is a quasinormal mode fluctuation of the black brane geometry. 

The key observations are as follows : 

1. Even if the horizon fluctuates in presence of the non-equilibrium energy-momentum and charge current fluctua- 
tions in the dual state, i.e. the metric and gauge field quasinormal modes in the background, in the perturbation 
expansion, we need to apply the incoming boundary condition and regularity only at the radial location of the 
horizon at late time, which in our coordinates is always at r — P /r^^. 

2. The quasinormal modes always have a negative imaginary part in their dispersion relation, so the pole in the 
complex frequency plane of the background perturbation will always be in the lower half plane. 

The first point above makes sure that we can always write the non-equilibrium part of the solution as the incoming 
and outgoing homogeneous solutions plus a particular solution completely specified by the source at the horizon 
exactly as in the case of the hydrodynamic shear mode. The second point will imply that integration over the 
background frequency will produce a divergence at the horizon unless we put the coefficients of both the incoming 
and outgoing parts of the non-equilibrium part of the solution to zero. Therefore, the non-equilibrium part of the 
solution is completely determined by the equilibrium part of the solution for any background quasinormal mode. We 
can thus simply repeat the exercise as we have done for the hydrodynamic shear-mode perturbation to obtain the 
retarded propagator for any background quasinormal perturbation. 

One may wonder if our prescribed solution at the horizon involving the specific particular solution is itself regular 
at the horizon. We have checked this is always so for the scalar field. In case of the fermionic field, we can repeat the 
arguments we have made in case of the hydrodynamic shear-mode. 
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For instance, let us consider a quasinormal mode for metric perturbation in the tensor channel with momentum 
k(f,) = 0. The frequency will be complex with a negative imaginary part as in (^). The explicit metric and gauge 
field for such a spatially homogeneous perturbation is as in (^T|). We can check that our prescribed non-equilibrium 
solution for the scalar field dies down at the horizon due to the factors : 



multiplying the equilibrium incoming wave solution with n and m being positive integers 
The general dispersion relation for a quasi-normal mode may be written as : 

W(b)(k(b)) = WR(b)(k(b)) - ia;i(b)(k(b)), with a;i(b)(k(b)) > 0. 



(99) 



Also both a;R(b) and a;i(b) admit Taylor expansion in k(b) (and do not vanish when k(b) — 0). The bosonic retarded 
propagator will take the following form in such a background : 



Gfl,(w,k,x,i) = J J '^'^^^ 



O(o)(c.o,ko) 



J(°)(wo,ko) 

^/0W(^ 
J(o)(t^o,ko) 2^^l O(o)(t^o,ko) 



5(c^-c^o)'5'(k-ko) 



O(o)(wo,ko) 1 /^0(iH^o,ko,k(b))^2^j^_j^^ 



j(i)(^wo,ko,k(b) 
J(°)('^o,ko) 



k-kn 



{uj-ujQ-\ (wR(b) (k(b) ) - it^i(b) (k(b) ) 



- t^o + 5 (wR,(b) (k(b) ) - iwi(b) (k(b) ; 



gik(b, -Xg"^ (k(6) )-icJi(b) (k(b) ) j * 



(100) 



The non-equilibrium part of the source and expectation values of the dual operators j'^-* (w, a;(b), k, kj-i,-)) 
and O*^^^ (w, W(b), k, k(b)) can be determined from the non-equilibrium part of the solution. J^^^^ (w, k, k(b)) 
and 0(inw,k, 

k(b)) appearing in the retarded propagator above are the residues of j'^^ (w, w^^,), k, k(]3)) and 
O'^"'^^ (cli, CLi(b)5 k, k(b)) respectively in a;(b) at a;R(b) (k(f,) ) — iwi(b)(k(b))- These will be linear in the hydrodynamic 
fiuctuations Sui, ST, Sp and the non- hydrodynamic fluctuations STr^J^\ vq and v'f^^\ and will have a systematic 
expansion in k^^) 

Similarly, the fermionic non-equilbrium retarded propagator will take the general form : 



G_R(a;,k,x,<) 



i I duj() I dfikc] 



2?(")(l^o, ko)7-<^(c^ - c^o)^'(k - ko) 



- ^ ( (wo , ko ) i'-TIa (wo , ko , k(b) ) 5^ (: 



k-ko 



-TZb (wo, ko, k(b)) (wo, ko)7-<5' (k-ko 



'(b) 



!;(b) 



gik(b) -Xg-' ("R(b) )-i^Ub) (k(b) ) j t 



w - Wo - 5 (wR(b) (k(b) ) - iwi(b) (k(b) ) 

1 

w - Wo + i (wR(b) (k(,,) ) - iwi(b) (k(b) )^ 



(101) 



This can checked by expanding /i(s,a;(„)) in (^) in Though this expansion as noted before is dangerous for seeing manifest 

regularity of the metric, it does good job for analyzing the behavior of the scalar field in the perturbed background. 
The Taylor expansion is ^(h) always make sense near equilibrium as the perturbations are slowly varying in space. However, all time 
derivatives need to be summed up for non-hydrodynamic perturbations at each order in the amplitude and ^{h) as the variation of these 
modes in time is not small even near equilibrium. 
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TZa and TZ b can be determined from the non-equilibrium part of the sohition via the defining relations : 

0(^)(^a;,k,tj(b),k(b)) = 7^A (w, W(b), k, k(b)) 0^^\uj,k), J^^^ (^w, k, W(b), k(b)) = (w, cj(b), k, k(b)) J^''^(w,k). 

(102) 

7?.^(ti;, k, k(b)) and 7?.B(a;, k, k(b)) appearing in the retarded propagator above are the residues of 7?.a(w, wjb), k, k(b)) 
and 7?._B(w,W(b),k,k(b)) respectively in i:j(b) at WR(b)(k(b)) - iwi(b)(k(b)). Both 7?.a(w, k, k(b)) and 7^b(i:j, k, k(b)) wih 
be linear in the hydrodyanamic fluctuations (5u,;, (5T, 5p and the non-hydrodynamic fluctuations 5'k\^^\ vq and 
and will have a systematic expansion in k(b) ■ 

Thus we indeed obtain an universal form of the holographic non-equilibrium retarded propagator (and hence the 
spectral function) in linearized non-equilibrium backgrounds at sufficiently late time. 



IV. NON-EQUILIBRIUM FERMI SURFACE AND DISPERSION RELATIONS 

We will show here that our prescription for obtaining the non-equilibrium retarded correlator gets a lot of support 
from field theoretic comparisons. We will begin with a brief review of how we obtain non-equilibrium correlation 
functions in field theory. Then we will show how our prescription reproduces the strongly coupled version of non- 
equilibrium dynamics at the Fermi surface in Landau's Fermi-liquid theory, and the non-equilibrium modifications of 
quasi-particle dispersion relations expected in field theory. 



A. Comparison with field-theoretic approach 

In field theory, there is no partition function which can play the role of a generating functional of non-equilibrium 
correlation functions. The way we obtain these is to construct a generalized effective action T {Oi{x),Gw{x,y)) 
whose arguments are not only the expectation value of the operator but also the two-point correlation functions 
of the operators. Extremizing this leads us to obtain non-equilibrium correlation functions as functionals of the 
expectation values of the operators in equilibrium and non-equilibrium states. The crucial point is that the generalized 
effective action has no dependence on temperature or other equilibrium/non-equilibrium parameters It is defined 
as a double Legendre transform of a vacuum observable constructed over the Schwinger/Keldysh closed real time 
contour as briefly reviewed in appendix D. Both equilibrium (temperature and chemical potential dependent) and non- 
equilibrium dynamics of expectation values of operators and their correlation functions can be derived by extremizing 
this generalized effective action. At equilibrium, we can take an alternative route by constructing a generating 
functional of thermal correlation functions as in vacuum, but in order to obtain non-equilibrium correlation functions 
the use of the generalized effective action is indispensable. 

We would like to mention here that the generalized effective action not only allows us to obtain the non-equilibrium 
two-point correlation functions, but it is also sufficient to obtain the three, four and higher point correlation functions 
. This is possible because through the effective action, we know the two point correlation function as a functional 
of expectation values of operators, i.e. we know them not in one but in a manifold of states. Furthermore, the effective 
action technique readily ensures that we satisfy Ward identities. In practice, we need to make an uncontrolled but 
educated approximation of the effective action which allows us to obtain non-equilibrium dynamics of expectation 
values of operators and their correlation functions. This has been successful for instance in the case of dilute cold 
non-relativistic Bose gases in optical traps [^0| , and in constructing a quantum kinetic theory of hadrons for modeling 
their evolution after their chemical and thermal freeze-out in the RHIC fireball 

The important point to note is that we can obtain the non-equilibrium correlation function by extremizing the 
effective action with respect to the correlation function first as below : 

Thus we obtain the two point correlation functions as functionals of expectation values of the operators. Here the 
time contour is the Schwinger-Keldysh closed real time contour, so this determines both the statistical function and 



This is also true for kinetic equations, like the Boltzmann equation. These equations do not depend on temperature or non-equilibrium 
parameters, the latter parametrize equilibrium and non-equilibrium solutions of these equations. 
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the retarded propagator (or the spectral function). Further when we substitute the extremal forms of the two-point 
correlation functions in the generalized effective action, we obtain the ordinary effective action, i.e. 

T{Oi,G'}AOi))=T{Oi). (104) 

Extremizing this further we obtain non-equilibrium dynamics of expectation values of operators. 

It is certainly interesting to see if we can construct a generalized effective action to obtain non-equilibrium correlation 
functions in holography too. This will allow us to determine the statistical function also and not the retarded 
propagator alone as we have done here. Work in this direction will appear in pO| . However, we note two crucial 
points of our holographic prescription for obtaining the retarded correlator (cquivalently the spectral function). 

1. Our prescription obtains the non-equilibrium retarded propagator as a functional of the expectation value of 
the energy-momentum tensor and the charged current parametrized by T, p, ST, Su, dp, and Trj"''''. 

2. The non-equilibrium part of the correlation function is determined completely by the equilibrium part through 
universal rules at the horizon which do not depend on the non-equilibrium state concerned. The rule simply 
involves putting the homogeneous pieces of the non-equilibrium part of the solution of the bulk bosonic/fermionic 
field to zero at the horizon. 

Putting these together, we can see a parallel with field theory. In both approaches, we do not need a specific rule for 
each non-equilibrium state, there is a universal rule which allows us to extract the non-equilibrium correlation functions 
from observables defined at equilibrium. In field theory the equilibrium temperature arises as the boundary condition 
appearing in the far future. The generalized effective action as mentioned before is just the double Legendre transform 
of an equilibrium observable, therefore non-equilibrium dynamics can be obtained from equilibrium observables in 
field theory as well. Furthermore, our holographic prescription has the same measure of universality as the generalized 
effective action to bring all non-equilibrium spectral functions under one fold at least in perturbation theory. 

The advantage of the holographic approach is that the late time behavior of the non-equilibrium spectral function 
is reproduced automatically without any need for resummation. Thus we can do conventional perturbation theory. 

B. Non-equilibrium dynamics at the Fermi-surface 

It might have been a bit surprising that the logic of regularity required that we put the extra boundary condition 
needed to determine the non-equilibrium part of the solution completely, at the horizon instead of at the boundary. It 
might seem that it would have been more natural to suppose that the source does not fluctuate from it's equilibrium 
value, so a Dirichlet condition at the boundary would have been more justifled. As we have already argued, this is 
not the case - the source gets screened or dressed by the collective excitations present in the non-equilibrium state 
also. From the holographic perspective, the horizon determines the screening/dressing of the source. 

We will here give another holographic interpretation of the non-equilibrium modification of the source. This will 
further vindicate that we need to put the extra universal boundary condition at the horizon and not at the boundary. 
That we have allowed the source to fluctuate from it's equilibrium value, is what will bring out the non-equilibrium 
oscillation of the energy per particle at the Fermi surface and non-equilibrium shifts in the quasi-particle dispersion 
relations. 

A hallmark of Landau Fermi-liquid theory is that the collective modes as captured by the Boltzmann equation leads 
to non-equilibrium dynamics at the Fermi surface. This dynamics is characterized by shifts in energy per quasi-particle 
at the Fermi surface Se at a given direction n and at a given point in space and time in response to a local fluctuation 
in occupation numbers of quasi-particles at the Fermi surface Sn. Landau postulated the following phenomenological 
relation @: 

5e{kph, X, t) = e{kph, x, t) — eo{kFh) = ^ /(n, n') Sn{kFn', x, t), (105) 

where eo{kFh) is the equilibrium energy of a quasi-particle at the Fermi surface which is just kp/2m*{T, ^) with 
m * (T, fi) being the effective mass at the Fermi surface dependent on temperature and chemical potential. The 
parameters /(n, n') are phenomenological inputs of the Landau model which can be obtained from field-theoretic 
two-point density correlation functions. These phenomenological parameters determine all thermodynamic and many 
non-equilibrium properties of Fermi liquids. 

To obtain non-equilibrium properties one has to assume validity of Boltzmann equation for Sn. The equilibrium 
distribution n^^^ is the Fermi-Dirac distribution at a fixed temperature and chemical potential and is a trivial solution 
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of the Boltzmann equation. Using (105) and the Boltzmann equation, it can be shown that the fluctuations 5n follows 



/(n,n) 



= / 



d6n{kpn,x,t) kph dSn{kpn,x, 

dt m*{T,^) dx 9k ^ dx 

(106) 

in the linearized limit. Above / captures the so-called quasi-particle collision kernel. Studying this equation we 
can extract all collective excitations including the zero sound, hydrodynamic shear-mode and non-hydrodynamic 
relaxation modes. In order to obtain the zero sound velocity, the collision kernel is n ot necessary but it is so in order 
to obtain the viscosity and relaxation modes. Substituting a solution for Sn in ( 105 ) we can obtain the oscillation of 
the energy per particle at the Fermi surface. 

The crucial point is that the oscillation is related locally to the fluctuation in the occupation number of the quasi- 
particles in (105). So, the oscillation in energy per particle at the Fermi surface is in sync with the propagating 
collective excitation. 

We note that in non-equilibrium, we cannot obtain the change in energy at the Fermi-surface by looking at the 
spectral function alone. This is because the non-equilibrium change in the spectral function comes from both (i) 
the shift of the residue, and (ii) the shift in the pole itself. We need to identify which part of the non-equilibrium 
contribution comes from the shift in the residue and which part comes from the shift in the pole. Moreover, the 
situation could be worse, as there can be non-equilibrium contributions which are simply analytic near the location 
of the equilibrium Fermi surface and be neither the shift in the residue nor shift of the pole. 

In the holographic set-up, the Fermi surface(s) is related to the existence of normalizable mode(s) of the bulk 
fermion field at zero frequency on a fixed momentum shell jist . As the black brane retains rotational symmetry, the 
Fermi surface is spherical (circular for a 2 + 1 dimensional system) . We will be working in 2 + 1 dimensional system 
(i.e. in a 3 -I- 1 dimensional bulk) for the sake of concreteness. 

It will be worthwhile for us to first define the Fermi surface holographically in a more general background which 
may not preserve rotational symmetry. This will help us to readily understand non-equilibrium dynamics at the Fermi 
surface. 

A Fermi surface picks up an internal direction in spin space. Therefore, let us represent first an arbitrary normalized 
complex 2-vector which picks up a direction in spin space by two real angles 9 and as below : 



cos6'e*'^ 



smf 



(107) 



The vector above may still be multiplied by an overall phase, but this will be unimportant for us. We then note that 
the hcrmitian matrix P defined as 



P(0,0) = 



.c2, 



cos tl sm ti e 



cos f sm U e 



-i20 



(108) 



is a matrix such that 



p2 ^ p p I cos tie 
smt/e" 



f COS tie 
I sin 9 e~ 



P 



smfe"' 
-cos9e~ 



= 0. 



(109) 



Therefore P is a projection operator, and it projects in the direction (107) and in the orthogonal direction it has 
eigenvalue zero. 

The holographic definition of Fermi surface at equilibrium is as follows. Let us choose a direction n in momentum 
space. Then there exists 9, (j) specifying a vector in spin space and kp for every n such that : 



P[9,(j)),GR{Lo = 0,k = fc^n 
P{9,(j))j(u} = 0,k = fc^n 



- 0, 
= 0. 



(110) 



where P is as defined in (108)and J is the source obtained from the bulk solution. The first equation above says that 
Gr is diagonal in spin space in the following basis : 



/ cos 9 e"^ 
\sin9e~^'t' 



sm t/ 

- cos 9 



(111) 
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which is the same basis in which P is diagonal. Thus this defines 9 and 4>- We note if we replace 9 hy 9 + tt/2, we 
merely exchange the eigenbasis. Therefore, if 6* is a solution, so is 6* + 7r/2. The second equation is equivalent to : 



J[lj = 0,k= kpiij ^ i^^'^j , cos9e-"^xi ^ smee"^X2- (112) 

Thus we have one linear complex equation to define kp. Therefore kp is complex (at finite temperature) and associated 
with a specific direction in spin space. To get the Fermi surface associated with the orthogonal direction in spin space 
which is also an eigenvector of P and Gu we need to solve above with 9 replaced hy 9 + 7r/2, i.e. 

j(uj = 0, k = kpnj = (j^^^ , sin6'e-*"^xi = - cos9e"*'x2- (113) 

As the AdSi Reissner-Nordstorm black brane background preserves rotational invariance, 9, <j) and kp will be inde- 
pendent of n. 



More generally, the holographic Fermi surface is kp{n) which solves ( 110 ) and is associated with a specific direction 
in spin space in which the retarded propagator can be diagonalized. The general definition stated here should be 
useful in analyzing cases where we have spontaneous symmetry breaking in the boundary, particularly when these 
order parameters break rotational invariance [Q. We note that at zero temperature kp is strictly real and corresponds 
to the pole at a; = 0, but for non-zero temperatures the pole at a; = is complex. The imaginary part of the pole 
is negative and represents smearing of the Fermi surface at finite temperature, and vanishes as the temperature is 
reduced to zero. Thus we can think of fc^ as a complex parameter whose imaginary part vanishes at zero temperature 
and has a small T expansion. The real part oi kp also has a small T expansion and is the Fermi surface. There is no 
dependence on uj as to find the Fermi surface uj is set to zero. In the Reissner-Nordstorm black brane, the dependence 
of the negative imaginary part of this complex parameter kp on the temperature is given by a power law for small 
temperatures 43 . This power is controlled by the near horizon ^^5*2 x geometry. 



It can also be shown that the retarded propagator and the spectral function also have a pole precisely when the 
source vanishes. Therefore, the holographic definition of the Fermi surface matches with the conventional definition 
which is that it is the location of pole of the spectral function in momentum space at vanishing frequency. In 
holographic systems we typically get a family of nested Fermi surfaces. 

As an aside let us mention that the pole structure of the holographic spectral function at equilibrium is different at 
small frequencies from that of a conventional Fermi liquid and the scaling exponents are controlled by the near-horizon 
AdS2 X i?^ geometry [|l^. This means that holographic systems have generically non- Fermi liquid behavior. 

The full non-equilibrium source is : 

j[^,t)=Jd^x (^J(")(c.,k) + j(i)(c.,k,k(b))e*(w-e"^("^'^'<'''^''-'"^<^'^'^^^ (114) 

We recall that the full source J can be determined from the boundary behavior of our prescribed solution for ^_|- 
through (|87|). In fact we can explicitly write in case of the hydrodynamic shear- mode up to first order in the 
hydrodynamic momentum k(h) : 

jW(c^,k,k(h)) = Ji'^(c^,k)<5u(k(h)) - k + 4'^(c^,k) fc,fc,fc(h),fc,(k(h)), (115) 

where J^'^ and Jg^ can be determined from the solution. 

We will be interested in obtaining the energy oscillation Suj{h,x,t) at the Fermi surface by calculating shift of the 
frequency pole for a fixed Fermi momentum. We have to solve this perturbatively in the momentum of the collective 
background mode k^i,) . 

Perturbatively, the energy shift on the Fermi surface Suj in the direction n at a given point in space-time is thus 
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obtained by solving : 

6uj{n, X, t) ^P(6l(°\ (a; = 0, k = kpnj ^ + 60{kFn, x, t) ^9eP(6l(°\ 0*°^) (cj 0, k = kpn 

+50(n, X, t) |^a^P(0(o) , </.(o) ) J(0) = 0, k = fc^^n) ^ = - |^P(0(o) , (c^ = 0, k = fc^^fi, k(b) ) ^ 



ik(b)-x 



~« (^t^R(b) (k(b) )-i;^i(b) (k(b) )^ t 



e V / , (116) 



where 9^-^^ and cj)^^^ label the spin orientation of the equilibrium Fermi surface as discussed before and P is as defined 
in ( |108| ). The above amounts to two complex equations and we have four unknowns, namely real 69 and 5cj) giving 
change in the orientation in spin space and complex Su). As we have mentioned earlier, the change in orientation in 
spin space cannot be directly read off from the change in retarded correlator due to the ambiguity in identifying which 
change is due to shift in the pole and which change is due to shift in the residue. We can obtain the non-equilibrium 
shift in spin space at the Fermi surface from the non-equilibrium source directly. 

The shift in the energy of the equilibrium Fermi surface with orthogonal spin orientation can be obtained by solving 
the above equation with 9^-^'^ replaced by 9'^'^^ + Tr/2. 

Clearly in the hydrodynamic shear wave background, Suj takes the form : 

(5tj(n,x,t) = ^(5a;^(^n, fcF)'5u(k(h)) • n -|- Jws^n, fc^^ njnjfc(h)j'5Mj(k(h))^ e'''('''''^e~'4¥r (117) 
Therefore, we find that the holographic Fermi surface indeed oscillates in space and time in sync with the background 



collective excitation. Nevertheless in order to obtain the analogue of (105) in holography linking the spectral shift at 
the Fermi surface to the statistical shift (i.e. shift in occupation number) we need to obtain the statistical function 
holographically also. We leave this for the future. 

C. Non-equilibrium shifts in energy and spin of quasi-particles 

Not only the energy per particle at the Fermi surface but other normalizable modes with non-zero frequencies also 
receive space-time dependent shifts in energy at a given momentum in sync with the background collective excitation. 
This can be interpreted as the space-time dependent shifts of the dispersion relations of the quasi-particles in the 
non-equilibrium medium. This is certainly expected as quasi-particles receive a thermal mass and if the temperature 
oscillates for instance, the dispersion relations indeed become space-time dependent. This is usually a hard calculation 
in non-equilibrium quantum field theory, but we can readily generalize the holographic strategy discussed above to 
obtain non-equilibrium shifts in quasi-particle dispersion relations. 

A particular quasi-particle branch can be identified via the following steps at equilibrium. 

1. Consider the equilibrium Green's function G^''(a;,k). This can be diagonalized at a given uj and k and the 



eigenvectors can be labelled as in ( |TTl| ) by 9'-°\uj,'k) and (f>^°\uj,k). Furthermore, if 6'(°) is a solution, so is the 
9^^'^ + n/2 as this merely exchanges the eigenbasis. 

2. The quasiparticle pole a;*-°^(k) can be identified with a definite orientation in spin space by solving : 

P(0(°),0(°)) J(")(c^(")(k),k) =0. (118) 

The above amounts to one complex equation which determines w^*'-'(k) with ^'■"^ and 0^*^^ determined in the 
previous step. The imaginary part of w'^'^^(k) is negative. To obtain the quasi-particle branch with opposite spin 
orientation, we need to solve the above with 9'^^^ replaced by 9^^"^ + 7r/2. 

Once again, if there is rotational symmetry in the background, i.e. if there are no order paramaters of spontaneous 
symmetry breaking which breaks rotational invariance, 9^^\ 0^"-' and aj^°''(k) can depend only on the modulus of k. 
The space-time dependent shift in dispersion relation is characterized by : 



u ^ UJ 



(k)+fo(k,k(b),x,i). (119) 
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The shift 5u: can be obtained by solving : 

(5c.(k, X, t) |^P(0(°) , )9<.J(") (u = (k), k) ^ + 6e{\,, X, t) {dgPid^''^ , J(°) (u = J"^ (k), ] 

+<5</)(k, X, t) |^a^P(0(o) , </)(") )J(") (c. = (k) , k) j = - |^F(0(") , ^(0) ) J(i) (c. = (k), k, k(b) ) ^ 

g«k(„.x ^-«(-'R(b)(k(„)-»c.:(„(k(„))t^ ^^^^^ 

The above equation amounts to two complex equations which allows us to solve the real unknowns 69 and S4> giving 
shifts in spin space and the complex unknown 6u). To obtain the non-equilibrium shift in the dispersion relation 
for the other equilibrium branch with orthogonal spin orientation, we need to solve the above with replaced by 

The solution of duj will take the form in a hydrodynamic shear-wave background, for instance, clearly takes the 
form : 



(5w(k,x,t) = ^(5wA(k)(5u(k(h)) • k + ^^^(k) fc,fcjfc(h),^Uj(k(h)) j e*''<")-^e-'^*. (121) 

Therefore, we see that the shift in the dispersion relation of the quasi-particle pole is also in sync with the propagating 
collective mode. Furthermore, though we have discussed the fermionic case explicitly here, clearly the same strategy 
can be applied to the bosonic field also. In fact, the source being a complex number instead of a complex 2-vector in 
the bosonic case, the equations will be much simpler. 

The shift Suj in the quasi-particle pole is generically complex. Interestingly the sign of the imaginary part of Suj can 
be both positive and negative. Thus we can get both non-equilibrium suppression or enhancement of quasi-particle 
decays as indeed observed in the RHIC fireball for various resonances (short-lived quasi-particles) |10[. 
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V. TAKING INTO ACCOUNT NON-LINEARITIES IN THE DYNAMICS OF THE 

NON-EQUILIBRIUM VARIABLES 

It is known that solutions of gravity which have regular future horizons reproduce non-linear phenomenological 
equations for irreversible processes in the dual field theory. The best studied examples are related to fluid/gravity 
correspondence. The full non-linear Navier-Stokes' equation with higher derivative corrections can be reproduced 
from gravity and this success has also been extended to the case of charged hydrodynamics j|] . As we have discussed 
before, gravity is expected to reproduce the general phenomenological equations which describe the full evolution of 
energy- momentum tensor and conserved currents which generalize hydrodynamics Q . This has been shown explicitly 
for the case of spatially homogeneous relaxation [Q. In all cases, the regularity of the future horizon determines the 
phenomenological coeHicients. 

We would like to show that the prescriptions described so far can be readily generalized to include non-linearities 
in the dynamics of the energy-momentum tensor and conserved currents characterizing the non-equilibrium states. 
We can systematically include these non-linearities into the retarded correlator, the shifts in the dispersion relations 
of quasi-particles, etc. 

The key is to see how the solutions for the bosonic and fermionic fields get determined in the perturbed background. 
Let us focus on the case of the hydrodynamic background. If we take into account non-linearities in (5u(k(h-)) in the 
background, clearly these non-linearities will also appear in the Laplacian of the bosonic field. Let us consider 
quadratic dependence on two distinct velocity perturbations (5u(A;(h)) and 6u[k'i^^^) for instance, at a given order in 
the derivative expansion m (i.e. at the mth order in the hydrodynamic momentum). The solution for $ will receive 
a correction quadratic in the amplitude of velocity perturbation and at mth order in the derivative expansion which 
can be represented as: 

$(2'™) (r, k, fc(h), /fc(h)) e'C^+'^f) )-^. (122) 
The radial dependence above can be determined the equation of motion : 

ai^^S^k - k')<i>^''^-\r, k', fc(h), fc[i,)) = fc(h), (123) 
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where D^^^ is the Laplacian for a scalar with 3-niomentum k in the unperturbed ^^5*4 Reissner-Nordstorm back- 
ground and S'*^^'™' is a generic source term. For m = 1 the source S*^^'^^ can contain terms hke (k • (5u(A:(ii)))(k(h) • 
(5u(fc^jjj)), etc. It also contains the solutions at the lower order in the perturbation expansion for instance 

Clearly the general solution of $ near the horizon can again be separated into two homogeneous pieces, the incoming 
and the outgoing modes, and a particular piece which has no arbitrary integration constant and is completely deter- 
mined by the source term 5^^'"*^. In order to satisfy the incoming boundary condition, we should put the coefficient of 
the outgoing mode to zero. Also as discussed before, the integration over the hydrodynamic frequencies a;(h) and w^j^^ 
will produce a divergence at the horizon for the incoming mode, as for instance in the case above with dependence 
on two hydrodynamic shear wave background modes like : 



rroV 
/2 J 



(h) '"(h) 



(124) 



Obviously the coefficient of the incoming mode has to depend on 6u and the hydrodynamic momenta required by 
the order in the perturbation expansion. The contribution from the frequency pole in 6u{uj(yj,'k(^^) given by the 
hydrodynamic shear dispersion relation produces the above divergent behavior. In general the divergence will always 
be there for any quasinormal wave background as it's dispersion relation a;(h)(k(h)) will have a negative imaginary 
part. Therefore, we should put the coefficients of the incoming mode at the horizon to zero too. Wc are just left 
with the particular piece which is completely determined by the source term containing the solutions at the lower 
orders in the perturbation expansion. Therefore, applying induction, at each order in the perturbation expansion, the 
non- equilibrium solution is uniquely determined by the equilibrium solution, i.e. the solution at the zeroth order in 
the unperturbed black brane background. The consistency of holographic duality requires the solution at each order in 
the perturbation to be regular at the horizon. 

As the solution is uniquely fixed at each order in the perturbation expansion, we can obtain the non-equilibrium 
contributions to the source and the expectation value of the dual operator by studying the asymptotic behavior of 
the solution at each order. This procedure can also be applied for fermionic fields. 

Once the source is obtained at a given order in the perturbation expansion, it is straightforward to obtain the shift 
in the dispersion relation of quasi- particles. For example, Su)^^'^'^\k.,x.,t) along with the non-equilibrium shift in the 
spin orientation given by 56^^'"^\\i,x,t) and 6(j)^'^''^\\i,x,t) can be obtained from by solving : 

X, t) (^P{0<-°^ , <^(°) )a„ J(°) = (k), k) j + 5^(2,™) X, t) (^d0P{0<-°^ , (i><-°^)J^°^ = (k), k) j 
+6<j)^^'"'\k,x,t)(d^P{0^°\^^°^)J^°^(^uj=u^°\k),kU = -(p(0(o),0(o))j(2.™)(a;=w(o)(k),k,k(b)) Je^(''(-)+'''(-))-^ 



( (i^RCb) (k(b) )-l-WR(b) (k' (b) )) -»(wi(b) (k(b) )-l-wi(b) (k' (b) ) j « 



e ^ ' . (125) 

A consistent perturbation theory for the solution in the non-equilibrium background thus suffices to take into 

(n/ 



account non-linearities in (5uj, (5T, ^I'Klf'^ etc. in the retarded correlation function, spectral function, non-equilibrium 



shift in the dispersion relations, etc. 



VI. DISCUSSION 



This paper has been devoted to developing a general holographic formalism for determining non-equilibrium retarded 
correlator, spectral function, shifts in dispersion relations, etc. Needless to say, we would like to use this formalism 
to numerically calculate these space-time dependent quantities in the specific set up of charged bosonic and fermionic 
fields minimally coupled to Einstein-Maxwell gravity in AdS\ discussed here. In particular, the following questions 
require attention. 

1. It is known that at equilibrium the temperature modifies the spectral function only in the infrared, while in the 
ultraviolet the spectral function remains as in the vacuum. It can be expected that we have a similar feature 
even in non-equilibrium - the ultraviolet behavior of the spectral function, quasi-particle dispersion relations 
should be independent of the state. It will be interesting to see if this is reproduced in our case. Some of 
the background quasi- normal modes indeed can have very high frequencies, while high frequency dependent 
corrections can also be generated by non-linearities. Therefore, numerical studies can help us understand how 
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the effect of fiigh frequency dependent background modes gets suppressed in the ultraviolet, if this is indeed the 
case. 

2. The non-equilibrium shifts in the dispersion relations can have both positive and negative imaginary parts. If 
positive it leads to suppression and if negative it leads to enhancement of the decay. It will be interesting to 
see if one can use non-linearities to design a background in which a specific quasi-particle can be stabilized 
against decay to a large extent in a certain range of energies. This can allow us to observe otherwise short-lived 
quasi-particles. In particular, it will be interesting to see if some bound states of heavy quarks can indeed exist 
in the quark- gluon plasma at temperature 175 MeV. 

3. The quasi-particle dispersion relations can change non-analytically with the temperature particularly if there 
is level crossing. It will be interesting to design a non-equilibrium background where the temperature varies in 
space and time over the range in which this non-analyticity can occur and study exactly how the quasi-particles 
behave in such backgrounds. It will be interesting to learn from such holographic examples how to describe such 
non-equilibrium states in field theory. 

Work is in progress by the authors of this paper to tackle such issues numerically . Our prescription here gives an 
algorithm to tackle such questions in specific holographic models. 

Another direction we want to pursue in the future is to study non-equilibrium spectral functions in states corre- 
sponding to a plasma undergoing boost-invariant hydrodynamic expansion as in the RHIC fireball. This will give 
us insights into how hadrons are produced and transported in the medium, and finally get frozen chemically and 
thermally. 

Acknowledgments 

The authors would like to thank C. Gowdigere for collaboration in the initial stages of this work. The authors 
would also like to thank R. Gopakumar and G. Policastro for valuable comments on the manuscript. AM would 
like to thank C. Bachas and G. Policastro for inspiration and useful discussions. AM would also like to thank all 
the members of the AdS/CMT discussion group in Paris for inspiration. The discussion in section II. B has been 
motivated by questions posed by D. Mateos during AM's seminar at Department de Fisica Fonamental, Universitat 
de Barcelona. The results of section IV. B and IV. C followed from valuable discussions with G. Policastro. AM would 
also like to thank the Institute for Theoretical Physics, Utrecht University, the organizers of SCGSC 2011, London 
and the organizers of Rencontres Theoriciennes held at Institut Henri Poincare, Paris for opportunities to present this 
work prior to publication. SB would like to thank the organizers of Asian String School 2012 for their hospitality, 
and Theory Center, KEK, Japan and Yukawa Intitute for Theoretical Physics, Japan for opportunities to present this 
work prior to publication. SB would like to thank S. Mukherji and M. Natsuume for inspiration and useful discussions 
which have helped us improve our manuscript. The research of AM is supported by the grant number ANR-07-CEXC- 
006 of L'Agence Nationale de La Recherche. RI would like to thank the Division of Biology at Calfornia Institute of 
Technology, where part of this work was carried out, for their hospitality. 



38 



Appendices 

Appendix A: Eddington-Finkelstein vs Schwarzchild coordinates 

In order to see regularity at the horizon manifestly in the metric (|37| ) corresponding to hydrodynamic shear- 
mode perturbation of the AdS^ Reissner-Nordstorm black brane, we can consider the following change of coordinates 
following M : 



t ^ v + 



-*^(h)^fci (^) Su\k^y,^)e^^'^^-y^-"i-^^^ + 0{e'), (Al) 



where 



and 



k{a) = / d6^, (A2) 



h{a)^£db (i^^)fc(6). (A3) 



These new coordinates r, v and are ingoing Eddington-Finkelstein coordinates. 
In these coordinates, the metric assumes the form : 



dv dx" 



ro 



dx-' 



+0{e^). (A4) 
The bulk gauge field however no longer remains in the radial gauge and takes the form : 
. 1 Vigpro 



A. = %^(l-91)+0(."), 



A. = -:iSpi(i_!^)j„,(t,,,,),*,.,. -"„,.) (i_ii?i2ft(!^))+0(.')^ (A5) 

It can be checked that the gauge field is also regular at the horizon. A^, Ai vanish while Ar is a constant at the 
horizon. We can bring the gauge field back to radial gauge by a regular gauge transformation. 
Most importantly, the ij components of the metric is regular as 

— /i(a) — k{a) = terms which are regular at the horizon (i.e. at a = 1). (A6) 



So, the metric is manifestly regular up to the first order in the derivative expansion in these coordinates. 

We can implement this change of coordinates order by order in the derivative expansion. Even beyond the 
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fluid/gravity correspondence, such coordinate transformations can be implemented perturbatively to see manifest 
regularity |§. 

Appendix B: The general solution for the non-equilibrium profile of the scalar field 

At the zeroth order, the equilibrium solution for a given mode can be written as an arbitrary linear superposition 
of two linearly independent homogeneous solutions ^"^{k,r) and r). Here k denotes (i^,k) collectively. Thus 

$(°)(A:,r) = A(°)(fc)$^(fc,r) +B(°)(fc)$^(fc,r), (Bl) 

where A'^^'>{k) and i?'^°)(fc) are arbitrary. 

Using the method of variation of parameters, we can write the general solution for the equation of motion for 
the non-equilibrium part can be found and is as below : 

rr <&S(fc + fc(h),r')(vi+F2)(fc,fc(h),r')<f(°nfc,0 
<l>(fc,fc(h),r) = -<I)^(fc + fc(h),r-) I dr' ^ ^ 



-^^{k + k(^^,r) dr' ^ 1 (B2) 

% W[^^{k + fc(h),r'), $^(fc + fc(h),r')]r'2/(T?^) 



Above fc(ij) denotes (a;(ij), k(h)) collectively, W denotes the Wronskian of the two homogeneous solutions, and li and 
I2 are arbitrary setting the range of the two integrals. 

One can readily verify that the above is independent of the choice of and for fixed /i and Z2. To see the 



general behavior at the horizon given by ( 53 ) one can set to be $™ and $^ to be . 

Furthermore, one notes that the above is consistent with the derivative expansion for any li and I2 as the dependence 
on 5ui and fc(h) comes from Vi and V2 directly. Comparing ( |63| ) with ( p9| ) one gets that the explicit contribution to 

O^^^ and J^l^ comes from Vi, and the contribution to O^g and jjj"^' comes from V2. 

Appendix C: Vielbeins and spin connections in the hydrodynamically perturbed black-brane metric 



We calculate here vielbeins, their inverses (or einbeins) and spin connections for the metric (37) which corresponds 
to a black brane perturbed by a hydrodynamic shear mode. The notation we use here is the same as defined in section 
III.B of the paper. As noted there, to ease computations we will choose, without loss of generality that is in the 
y direction. As 5\i is transvere in the shear-mode, k(h) will be then in the x direction. On the other hand k can have 
arbitrary x and y components in order to retain full generality. 

The non-zero vielbeins upto first order of derivative expansion are : 



^, e^^- / . (Cl) 



From this, one can also construct inverse vielbeins (or einbeins) which are as follows : 
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e* = 



ei = - 



(5uj/(fc(h))e 



i('i:(h)xa:-W(h)t) 



e? = — 



1 f I rro 

21 



ef. 



■ 



(C2) 



In order to derive the equation of motion of the Fermions in the given background, we require the spin connections 

A B 

associated with the first order metric. The non-zero components of the spin connection, UJJ^~ are as below : 





3r3( 


-1 + / 




) 




^(h) 




//(^) 



tv 



yt 



ty _ 

OJr — 



tr 



yt 



rt 



4/2 



2P 



y y 

xy yx 





2Fi 



-izK(fc(h))e^^'="'"'^-""''*' (-l + /(^)) fc(h)x, 



Wjr- = — . 



6rro 
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yr 



yr 



ry 



ry 



(5%(fc(h))e'(^('>)-^-"('>)*) 


2/2(-l + /(l^))/( 


2^)-rro(l + /(^))/'(^)^ 


4Pr^ 







UJx — —Ulx — 



6rrQ 



U)y —I^V ~ 



/(^) 



(C3) 



Here prime denotes derivative with respect to 

It can be checked that the above spin connections satisfy Cartan structure equations up to first order in the 
derivative expansion. 



Appendix D: The generalized effective action 

We will review the formalism for bosonic operators here. The generalization to fermionic operators is straightfor- 
ward. 

The starting point of the construction of the generalized effective action is to generalize the partition function which 
is a generating functional of the vaccum correlation functions. Here on top of a source J; (a;) for a single operator 
Oi{x), we add a non-local source Kiii{x,y) for a pair of operators Oi{x) and Oi'{y), and define: 



Z{Ji,Kn') = e^^(^"^"') 



= J P$,exp J d''xJi{x)Oi{x) + ^ J d'^xd^'y Oi{x)Ku\x,y)Oi,{y) 

Above D is the number of space-time dimensions in field theory. 

We then define the expectation value of the operator Oi{x) and the Green's function Gii'{x,y) through : 



(Dl) 



6WiJi,Kii,) 
SJi{x) 

6W{Ji,Kw) 
6Ku'{x,y) 



Oi{x), 



-lOi{x)Oi,iy) + Gw{x,y)\. 



(D2) 



Eliminating J; and Kui in favor of O; and Gw , we can now do a Legendre transform to define the generalized effective 
action : 



T{Oi,Gu,) = W{JuKu-)- j d''xJi{x)Oi{x)~^- j d'' xd'^y Ku-{x,y){oi{x)Oi,{y) + Gn'{x,y)^. (D3) 



Clearly, 

6T{OuGu.) 



50i{x) 



^ = -Ji{x)- j d'^y Kir{x,y)Oi>{y), 



Therefore, in absence of sources, extremizing the generalized effective action V{Oi,Gii') gives the dynamics of both 
the operators and their Green's functions. 

Such an effective action is usually considered for the elementary fields and their Green's functions in the literature. 
However, as discussed above we can construct the same for the set of gauge-invariant single trace operators in a 
non-Abelian gauge theory. 
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There is one important point in the above construction. The effective action is constructed over the so-called 
Schwinger-Keldysh real time contour (Fig. which travels from — oo to oo infinitesimially above the real line and 
then back from oo to — oo infinitesimially below the real line. It is necessary to consider this "closed-time" contour 
because the usual time-ordered Green's function or the Feynmann propagator do not contain the full information 
about the operator in presence of sources in a non-equilibrium state as mentioned in the Introduction. The closed-time 
contour ensures we propagate the full information of the operator in presence of the sources. In fact, the full closed- 
time contour ordered Green's function can be written as a combination of the commutator and the anti-commutator. 
For instance, if both operators are bosonic the 

Gcw{x,y) = \{{Oi{x),Oi.{y)}) - '-{[Oi{x),Ov{y)]) ^igudx'' - y^). (D5) 

Above C denotes the closed-time contour, and x^ and y^ are the time coordinates of the D-dimensional position vector 
X and y respectively. 

c 

> 

^ 

< 

t 

FIG. 2: The closed time Schwinger-Keldysh contour is as above. The forward and backward directed parts of the contour have 
been displaced slightly above and below the real axis just to distinguish them clearly. 



In fact, as discussed in the Introduction, the spectral function is related to the commutator and the 

statistical function (or Keldysh propagator) Gjciv {x, y) is related to the anti-commutator in the following way (for 
bosonic fields): 

Aw{x,y) - i{[Oi{x),Oi,{y)]), 
GKw{x,y) = ^{{Oi{x),0,{y)}). (D6) 

The coupled equation of motion of the spectral and statistical functions are obtained from the generalized effective 
action. 

The generalized effective action has no dependence on temperature or non- equilibrium variables, it is defined as 
a Legendre transform of the vacuum persistence amplitude in the presence of single and double operator sources. 
However, thermal and non- equilibrium propagators also can be obtained as solutions which extremize this generalized 
effective action. In order to obtain thermal propagators, we need to impose translational invariance, so the Wigner 
transformed spectral and statistical functions Aw (w, p, x, t) and Gk:u' (i^, P, x, t) do not depend on the centre-of-mass 
coordinates x and t. Furthermore, they should be related by a temperature dependent fluctuation-dissipation relation. 
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